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^ ' Abstract 

<N : 

We study the Hamihonian system made of weakly coupled anharmonic oscihators 
arranged on a three dimensional lattice 'I-2N x and subjected to a stochastic forcing 
mimicking heat baths of temperatures Ti and T2 on the hyperplanes at and N . We 
I introduce a truncation of the Hopf equations describing the stationary state of the system 

. which leads to a nonlinear equation for the two-point stationary correlation functions. We 

prove that these equations have a unique solution which, for large, is approximately 
a local equlibrium state satisfying Fourier law that relates the heat current to a local 
' temperature gradient. The temperature exhibits a nonlinear profile. 
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1 Introduction 

Fourier's law states that a local temperature gradient is associated with a flux of heat J which 
; is proportional to the gradient: 

J[x) = -k{x)VT{x) (1.1) 

where the heat conductivity k{x) is a function of the temperature at a; : k{x) = k{T{x)). 

Fourier's law is experimentally observed in a variety of materials from gases to solids at low 
and at high temperatures. It also belongs to basic textbook material. However, a flrst principle 
derivation of the law is missing and, many would say, is not even on the horizon. 

The quantities T and J in (1.1) are macroscopic variables, statistical averages of the variables 
describing the microscopic dynamics of matter. A flrst principle derivation of (1.1) entails a 
deflnition of T and J in terms of the microscopic variables and a proof of the law in some 
appropriate limit. 

An example of an idealized physical situation would be a crystal occupying the region 
[0, N] X in R^. The crystal is heated at the two boundaries by uniform temperatures, Ti on 

^Partially supported by the Academy of Finland. 
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{0} X and T2 on {A^} x [R^. Then, for large A^, the relation (1.1) should hold, with corrections 
o{l/N) and VT = 0{l/N) = J. 

The proper description of the crystal would be quantum mechanical, but, being a macro- 
scopic law, (1.1) is expected to hold for a classical system as well and the quantum corrections 
are expected to be small except at low temperatures. An example of a classical toy model of 
such a system is given by coupled oscillators organized on a o?-dimensional lattice Z*^. 

Consider a subset A = [0, A^] x Z'^"^ of Z*^ and let a; G A index the dynamical variables, 
coordinates and momenta p^. The dynamics in the phase space IR^^ = {{qx,Px) | a; e A} is 
defined in terms of the Hamiltonian 

i.e. the Hamiltonian flow is given by the system 

Qx = Px (1-2) 
. - . (1.3) 

We can think of this model as describing atoms with unit mass with equilibrium positions at 
x e A and q^ being the deviation of the position of the atom indexed by x from its equilibrium 
{qx should of course be in R'^ but we simplify and take q^ G IR), while Px is the momentum 
(velocity) of the atom. 

The potential V should describe the forces between the atoms and possible restoring or 
pinning forces pulling q^ to the equilibrium q^ — 0. Let, for unit lattice vectors 

Vaqix) = g(x + e„) -q{x) 

denote the discrete derivative. Then, the interactions are described by a potential t/(Vg) 
whereas the pinning is described by W{q), and V — U + W. In the simplest case, U and W 
are local: 

W{q) = E^fe)' 
xeA. 

u(Vq) = E «(Vag(^)), 

x,x+ea&A. 

and, to model small oscillations, w and u are given by a low order polynomial: 

w{q) = aq^ + bq^ + cq^, 

and u similarity. 

In addition to the Hamiltonian dynamics of the system, we want to model the heating of 
the system at the boundary. The simplest way to model this is to add stochastic forces to (1.3) 
for X G 9A, see Section 2 for details. Then the deterministic flow (1.2, 1.3) is replaced by a 
Markov process {q{t),p{t)) and the Fourier law will be a statement about the stationary state 
of this Markov process. 

The first question one would like to answer is the existence and uniqueness of the stationary 
state. When Ti = T2 = T such a state exists and is a Gibbs state, given formally by 

^e-m<i,p)dq dp (1.4) 
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(see Sect 2), with {3 = 1/T. This is customarily refered to as the equilibrium state in contrast to 
the non-equihbrium situation Ti ^ T2. In the latter case, there is no simple formula like (1.4) 
and, indeed, in our setup, even the existence of a stationary state is an open problem. In the 
d = 1 case of a finite chain of + 1 oscillators the existence is proved provided the interaction 
potential U dominates the pinning one W (see [9, 10, 11, 12, 17, 13, 24, 25, 26]). In this case, 
uniqueness is also proved, i.e. the Markov process converges to this state as t ^ 00. 

Supposing that wc have a stationary state ^, let us formulate the statement (1.1). The 
Hamiltonian flow (1.2, 1.3) preserves the total energy H. In particular, if we write H as a sum 
of local terms, each one pertaining to a single oscillator: 

H = ^ Hx, 

then, under the flow (1.2, 1.3), 

a 

where the microscopic heat current j{x) will depend on py and Qy for y near x (see (7.2) for a 
concrete expression). Let also t{x) = ^pl be the kinetic energy of the oscillator indexed by x. 
Then, the macroscopic temperature and heat current in eq. (1.1) are defined by 

T{x) = t{x) 
J{x) = j(a^) 

where E^ denotes expectation in the stationary state. 

We do not attempt here to give a comprehensive review of the status of (1.1), but refer the 
reader to the reviews [6], [20] and [29]. There is also substantial amount of work on Fourier's 
law for fully stochastic models (i.e. where there is noise in the bulk too), going back to [16], 
[18], see eg. [7], [3], [4], [5], [14], [15]. 

The only rigorous results in our model are for the harmonic case where U and W are 
quadratic [27, 28]. In that case, Fourier's law does not hold: the current j'ls 0{1) as N ^ 00 
whereas VT = except near the boundary. If the model has pinning, i.e. W ^ and U or W 
are not harmonic, the law seems to hold in simulations in all dimensions [1]. In the unpinned 
W = anharmonic case, conductivity seems anomalous in low dimensions: k in (1.1) depends 
on N as N"' in d — 1 and logarithmically in d — 2. It is a major challenge to explain the a 
which, numerically, seems to be in the interval [1/3,2/5] (see [21], [22], [23] for theories on a). 

One way to try to get hold of the stationary state fi is via its correlation functions. Denote 
Ux = {(lx,Px)^ and choose U,W even for simplicity. Then, the Hamiltonian vector field (1.2, 
1.3) is a sum of a linear and a cubic term in u. Therefore, the correlation functions 

Gn{xu X2, ■■; Xn) = -E^M^j (g) ... (g) Ux„ 

satisfy a linear set of equations 

where An,Vn,Cn are linear operators coming respectively from the linear and cubic terms of 
(1.2, 1.3), and the term coming from the noise via Ito's formula. See Section 3 for a detailed 
derivation of these equations. 
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Such equations for the correlation functions are known as the BBGKY hierarchy in particle 
systems or the Hopf equations in turbulence. Although linear, the system (1.5) is intractable 
due to the appearance of Gn+2 in the equation for Gn- 

In this paper we will consider the situation where the equilibrium Ti = T2 Gibbs measure is 
close to a Gaussian measure. This holds if the anharmonicity in u and w (the coefficients b and 
c) is weak and the harmonic part in w (i.e. the pinning) is large (i.e. the Gibbs measure is far 
from critical), as we will assume. In such cases, we expect that the non-equilibrium measure is 
also close to a Gaussian. In such a situation, one can attempt a closure of the Hopf equations, 
i.e. to express the higher order correlation functions Gn+2 in terms of 0^ with m < n, thereby 
obtaining a finite set of equations for Cm, m < n. 

We will introduce such a closure and solve the closed equations. The simplest closure would 
be to write, for n = 2, the equation 

and set = 0, thereby obtaining a closed quadratic equation for G2. It turns out that this 
is too simple: the solution will be qualitatively similar to the one of the harmonic case. Our 
closure is done to the G4-equation by setting the connected 6-point function to zero. This is 
an uncontrolled approximation that we do not know how to justify rigorously. An analogous 
approximation was studied in both classical and quantum systems in [29] and has been used in 
[19] in a model similar to ours, but in a translation invariant setting and in one-dimension; our 
model, in one- dimension, was further studied, theoretically and numerically, in [2]. 

Our motivation for studying the closure equations in detail is on the one hand in the 
interesting picture of the local equilibrium state that emerges and on the other hand in building 
approaches that go beyond this approximation. Traditionally one arrives to such a closure in 
an appropriate limit, the "kinetic limit", which in our case means taking the anharmonicity 
proportional to A^~^ and rescaling distances by A^, see [29]. One then arrives to a Boltzman 
equation and after a further limit [29] to the Fourier law. Most of the structure of the stationary 
state correlations disappear in these limits. 

In our case we arrive to approximate expressions of stationary correlation functions whithout 
any limits (admittedly with no control of the corrections!) and then can study how the Fourier 
law emerges from these expressions. Other interesting phenomena emerge too. In particular, 
one expects the presence of very long range spatial correlations in the stationary state even 
though the equilibrium state has exponential decay of correlations. Although we do not demon- 
strate the presence of these long range correlations, because we obtain only upper bounds on 
the decay rates, not lower ones, we show how to handle the technical problems caused by these 
slowly decaying correlations, and that could be useful in other contexts. Finally, we believe 
some of the methods developed in this paper could be of use in trying to prove the existence 
of the kinetic limit and Fourier's law therein. 

The outline of the paper (to which the reader can return later) is as follows: in the next 
section, we define our model and we derive the Hopf equation (or BBGKY hierarchy) in section 
3. In section 4, we explain the particular closure that we will study, leading to our final equations 
(see (4.18)- (4.20) below). Section 5 is devoted to several changes of variables: we first apply 
a Fourier transformation, and then we introduce variables that could be called slowly and 
fastly varying, namely the one of which the non translation invariant part of of the correlation 
functions depends (slow) and the one related to the transtation invariant part (fast). Next, we 
outline our arguments and state qualitatively our main result (section 6). To prove the latter, 
we first derive (in section 7) identities satisfied by our equations, which take the form of current 
conservation equations, consisting of an energy conservation law and a number conservation 
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law (the presence of the latter being, to some extent, a consequence of our approximations, 
i.e. of our closure). We also write down the stationary states in the translation invariant 
case. These do not reduce, for our closed equations, to the usual Gibbs states, but depend 
on two parameters, the temperature, as one would expect, but also a "chemical potential", 
corresponding to the number conservation law. These conservation laws are related to the 
presence of zero modes in the linearization of our equations, which are discussed also in section 
6. In fact, the current conservation equations coincide with the projection of the full equations 
on the zero modes. In section 8, we define precisely the spaces in which our equations are solved 
and we state our main result in a more technical form. The solution that we shall construct is 
the sum of a modified stationary state, with coefficients (temperature and chemical potential) 
slowly varying in space, and a perturbation. 

The main technical problem that we face is that the nonlinear terms in our equations 
involve colhsion kernels that are delta functions (or principal values) (see (5.12)- (5.14) below). 
Since we want to solve our equations by using a fixed point theorem, we need to show that 
the nonlinear terms belong to the space that we introduced, and, because of the presence of 
the delta functions, this is rather technical. Section 9 is devoted to solving those problems, 
but most of the proofs of that section arc given in Appendix B. Another problem is that the 
linear operator in our equations is not invertible, because of the zero modes. In section 10, 
we show that our linear operator can be inverted in the complement of the zero modes. This 
uses the fact that this operator is a sum of a multiplication operator and a convolution. To 
show invertibility, it is useful to know that the convolution operator is compact and this in turn 
follows from Holder regularity properties that are proven in section 9. Finally, in section 11, 
we prove our main result, which consist in using the result of section 10 to solve the equations 
in the complement of the zero modes with the solution to the current conservation equations, 
i.e. of the projection of the full equations onto the zero modes. The first equations lead to 
Fourier's law, namely an expression of the conserved currrents in terms of the parameters of 
the modified stationary state (temperature and chemical potential ) . And the second equations 
determine the spatial dependence of those parameters. 

2 Lattice dynamics with boundary noise 

Let us define the model we consider in more detail. Instead of working in the strip [0, A^] x Z^~^ 
it is convenient to double it to the cylinder V = 1-2n x Z^^^ where 'Z2N are the integers modulo 
2N. The noise is put on the "boundary" {0} x T^'^ U {A^} x J.'^'^. We consider the phase space 
(g,p) e [R^^ i.e. q = {qx)xev, Qx = Qx+y for y = (2A^, 0), and similarly for p^. The dynamics is 
given by the stochastic differential equations 



dqx = 



Pxdt 





where 



^ xeV ^ ^ xeV 



(2.2) 



Ix = 7(<^xiO + ^xxn) 



(2.3) 
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and the random variables ^x{t) are Brownian motions with covariance 

Eix^y = 4.i6xy{Ti6,^o + T24iiv)t. (2.4) 

The Hamiltonian (2.2) describes a system of coupled anharmonic oscillators with coupling 
matrix o;^: 

x,yeV 

Our analysis requires that the Fourier transform u!^{k) of o;^ is smooth and 

u;'^(k) = + p{k) 

with p{k) = 0{k'^) as /c — > 0, and > 2||p||oo- Moreover we will need some regularity 
properties that will be checked explicitely for 

a;' = (-A + mY, 

i.e. 

d 

cu(A;) = 2^(1 -COS A;a)+m^ (2.5) 

a=l 

see the proof of Proposition 9.3. 

As explained in the Introduction, the equations (2.1) describe a Hamiltonian dynamics 
subjected to stochastic heat baths on the "boundaries" of V, with temperature Ti on the 
hyperplane Xi = 0, and temperature T2 on the hyperplane Xi = N. This defines a Markov 
process {q{t),p{t)) in the phase space R^^ and we are interested in the stationary states for 
this process. 

If the temperatures are equal, Ti = T2 = T, an explicit stationary state is given by the 
Gibbs state at temperature T of the Hamiltonian H. This probability measure is given as a 
weak limit 



vt — lim exp 

M-+00 Zm 



A ^ 4 

xGVm 



Uridp, dq) (2.6) 



where Vm is defined by \xi\ < M, i = 2, . . . ,d, and fiT is the Gaussian measure with covariance 

EpxPy = T6xy, Ep^qy = (2.7) 

Eq^qy = Tuj-\x - y) (2.8) 

For small A, the Gibbs measure vt is nearly Gaussian, with (2.7) still true and small 0{X) 
corrections in (2.8). It is very well understood via cluster expansions. Physically, the fact that 
the Markov process reaches the stationary distribution vt means that the heat introduced at 
the boundary spreads inside the system, which reaches equilibrium at temperature T. 

When Ti 7^ T2, things are very different. Even the existence of a stationary state is not 
known rigorously (not even in finite volume, M < 00). However, physically, one expects a 
unique stationary state v to exist. In this paper we assume this and inquire about the properties 
of V. In particular one would like to understand how the heat from the boundary now spreads 
inside the system: what is its stationary temperature distribution and what sort of fiux of heat 
exists in it. 
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3 Hopf equations 



Let us introduce a more compact notation for the stochastic differential equation (2.1). Denote 
{q,pf = u, K{u) = -A(0, {Vu)^ = {^^IxPxf and 77 = (0,0^- Then (2.1) becomes 

du{t) = ((A - T)u + A{u))dt + dri{t) (3.1) 

where 



A 

Define the correlation functions 



1 




Gn{xi,...,Xn,t) = E U^,{t) ® ... ®'U^„(t) G R^"^ 

By Ito's formula, we get 

Gn = {An — Tn)Gn + A„G„_|_2 + CnGn-2, (3-2) 

where 

A„ = A (8) 1 (8) ... (8) 1 + ...1 1 (8) ... (g) A 
and r„ is defined similarly. Moreover, 



AnGn+2 = ^^xi ® ■■■ ® Mu)^^ (g) ... (8) 

i=l 

GnGn—2 ^ ] Cx^xj Gn—2 (^1 ; ■ ■ ■ ; ; ■ ■ ■ ; ^ j ■ ■ -Xn) 



where the arguments Xi, Xj are missing. This defines linear operators from R,2(n+2)y _^ K^nv 
and R2('^-2)^ ^ R2nv respectively. C equals one-half the time derivative of the covariance of 
rj, i.e. 



with 




C 



Cxy = '^l^xyiTiSx^Q + 725a;^Ar). (3.3) 

Suppose z/ is a stationary state of the process u{t). Then equation (3.2) leads to a linear 
set of equations for the stationary correlation functions 

Gn{,xi,...,Xn) = j ®l:=iUxi i^{du) (3.4) 

i^n ~ ^n)Gn + AnGn+2 + CnGn-2 — 0. (3.5) 
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The equations (3.5) have the drawback that they do not "close": to solve for G„, we need to 
know G„+2- 

For A small, the equilibrium Ti = T2 measure is close to Gaussian. When Ti ^ T2 we expect 
this to remain true; however, the measure will not satisfy Ep^Qy = 0. 

We look for a Gaussian approximation to equation (3.5) for small A by means of a closure, 
i.e. expressing the Gn in terms of G2- Since G„ 7^ only for n even, the first equation in (3.5) 
reads: 

{A2-r2)G2 + A2G4 + C^0 (3.6) 
The simplest closure would be to replace G4 in (3.6) by the Gaussian expression 

^G2{xi,Xj) <^G2{xk,xi) (3.7) 
p 

where the sum runs over the pairings of {1, 2, 3, 4}. Equations (3.6) and (3.7) lead to a nonlinear 
equation for G2- It turns out that the solution to this equation is qualitatively similar to the 
A = case, i.e. G2 does not exhibit a temperature profile nor a finite conductivity. The only 
effect of the nonlinearity is a renormalization of uj. We will therefore not discuss this closure 
any further. 

The next equation is 

(A4 - 14)^4 + A4G6 + C4G2 = 

Write 

G4^Y.G2®G2 + Gl 
p 

and 

Ge^^G20G20G2 + ^G2®Gl + Gl 

p p' 

where G^ and Gg are the connected correlation functions describing deviation from Gaussianity 
and the sums run over the usual partitions of indices. After some algebra, we may write the 
first two Hopf equations in the following form: 

{A2 - + 22)^2 + A2GI + C = (3.8) 

(^4 - r4 + ^4)G'i + h{G2) + A^Gl = 0, (3.9) 
where the operators E2 and S4 are 

MG2)G2 = A2 E <^2 ^2 (3.10) 

p 

^4(^2)^4 = Y.KG2®Gl 
p 

and A4 means the following: G2 <8) Gl belongs to (R^^)®^ ® (R^^)®^; is a sum of terms 

a; = E 

i<j<k 
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where A4 acts with A4 in the spaces i,j,k and as identity in the rest. A4 then has at least one 
of the indices i,j,k equal to either 1 or 2. 
Finally, 

P 

where A4 is similar to A4, but with A4 acting on all the three factors G2. Exphcitely, denote 
i = {a, x) e {1, 2} X R^, so that u = (ui), U(^i^^) = q^, U(2,x) = Px- Then, 

4 

&(G2)nj2i3M = 6 X! Kj2j3H n iG^)iaja + («i ^ib,b^ 2, 3, 4) 

jijsji a=2 

with 

4 

^hmzH = --^^012 n ^aal^xixa- (3.11) 
a=2 

Equation (3.9) may be solved for G4: 

Gl = -(A4 - r4 - E4)-^(6(G'2) + A4G'^), (3.12) 

provided that ^44 — 14 — E4 is invertible. Substitution of (3.12) in (3.8) yields a nonhnear 
equation for G2 with dependence on Gg: 

{A2-V2 + T.2)G2+M{G2,Gl)+C = Q (3.13) 

with 

N{G2, Gl) = -A2(A4 - r4 + S4)-^(6(G'2) + A4G'^) (3.14) 

If we set Gg = in (3.13) we get a closed equation for G2. However, before defining the closure 
equation to be studied, let us discuss in more detail the term 6(6^2) ■ 

4 Closure 

The leading term (in powers of A and 7) in (3.14) is 

M' = -k2Al%G), (4.1) 

where we write, as we shall do from now on, G for G2, since we shall only deal with G2. Let us 
write this more concretely. 

To define the inverse of A.i wc define the stationary state correlation functions G„ as limits 
e — of the stationary state where, in equation (2.1), we add a term ep^dt, and in (2.4) 
a term AetS^y, i.e. we put a noise and a friction of size e everywhere. Then, the matrix A 
becomes 
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and, letting 
we have, 

Thus, 



R{t) = e*^. 



(4.3) 



poo roo 

-^-1 = / e'^^dt = / R{tfUt. 
Jo Jo 



r^^K^^')^)...... n(^(^)G).„,„ + (3 permutations), 

j a— 2 

with j = (j2, ja, j4) and, then, with i = (^2,^3,^4), 



J 1 



+3 E _E r *A.„(fl(()A)^_. n (fl(t)G) . {R{t)G 



J Jl»2«3 



a=2 



+ (i ^i' 



RecaUing the expression (3.11), some calculation gives 

Kpi^,y) = 6A2^«2 rdtj:{Rit)G)jx,zr 

JO ^ 

■ 3Ru{t,x- z)(R{t)G)^^{y,z) + (R{t)G)^^{x,z)Rp2(t,z - y) 



+ tr 



(4.4) 



where tr means that both a, j3 and x, y are interchanged and where the translation invariance 
of R was used. Since dtR — AR we have, from (4.2), 



and, using dtR — RA, we get: 



where 



dtRia — R2a 

Rp2 = —dtiRGo)pi, 



Go 





1 



(4.5) 
(4.6) 

(4.7) 
(4.8) 



Inserting (4.6) in (4.4), and integrating by parts, (4.4) becomes: 

Af^pix, y) = Hapix, y) + 6A^5a2(^^i Qi^^ zfuj-^{z -y) + tr 

z 

where we denote G\\ — Q and 

Nap{x,y) = 18A2(5„2lim H dtY.[(R{t)G) ix^z)]^ 
R,2{t, x-z) {R{t)G)^^{y, z) + {R{t)G)^^{x, z) (R{t)Go) ^^{z - y)] + tr (4.9) 
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The closure equation that we will study is the replacement of the exact equation (3.13) by 

{A2-r2)G+J\f{G)+C^0 (4.10) 

i.e. we drop the terms E2, S4, Gg, as well as the second term in (4.8). Before discussing the 
motivation for this approximation, let us write (4.10) explicitely. Let 

C = ( p ) (4.11) 

so, e.g. H{x,y) — Eq^Py and H^{x,y) — H{y,x) . Then, (4.10) can be written as: 

H + H^^O (4.12) 

P -uj^Q - HT +Mu = (4.13) 

-uj'^H - H^uj^ - PT -TP + J\f22 + C = (4.14) 

where 

^xy = 7<^a;j/((^a;iO + ^xiAf) (4.15) 



and C was defined in (3.3). 
Let 



J^^(H-H^) (4.16) 



so by (4.12), H — J, — — J, and we write: 

G=(5j p)- (4.17) 
Then (4.13) and (4.14) can be written as: 

2P = u^Q + Quj^ + JT-rj- Hi2 - NI2 (4-18) 



'Q-gu;2 + rj + jr + A/i^-7Vi2 = o (4.19) 



u'^J - Jo;^ + PL + TP - jVas - C = (4.20) 

(4.18) - (4.20) is a system of nonlinear equations for the correlation functions J, P. 
An important property of TV" in equation (4.9) is that, for all T, 

U{TGq) = (4.21) 
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for 





and, for the equilibrium case, with 7 7^ 0, Ti = T2, only one of these persists, namely the one 
with T = Ti. 

Note that the true equilibrium Gibbs state has P — T, J — 0, and 



with (T = 0{X). The terms S2, S4 and the second term in (4.8) would contribute to changing the 
(T = of our closure solution to a o" which would agree with the true a to O(A^). Dropping these 
terms, as well as is done for convenience and should not change our analysis qualitalively. 
The uncontrolled approximation consists in dropping Gq. Presumably, for small A, this would 
not make a qualitative difference. It would be interesting to try to prove this in the kinetic 
hmit A = 0{l/y/N). In the rest of this paper, we will discuss only equations (4.18)-(4.20). 

5 Changes of coordinates 

For large we expect the solution to the equations (4.18)-(4.20) to be translation invariant 
in the directions perpendicular to the 1-direction with a slowly varying dependence of the first 
coordinate Xi. It is therefore convenient to represent G in coordinates that are suited to such 
behaviour. 

We first write eqs. (4.18)-(4.20) in terms of the Fourier transform of G. Recall that a; e Z'^ 
with xi G Z2N- Introduce momentum variables: 



Q{k) ^T(u;{k)^ + a{k,X,T)) 



-1 



q = (q,q), 



and write 




with the shorthand notation 




(5.1) 



where q e ^'Z,2n and / = [0, 27r]'^ ^. Then R becomes a Fourier multipher 
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where uj{q) = (<^^(g) — x)^^^- Then, letting e — >• in the matrix, we get 

R(t q) = ( cosa;(g)i ^ sin uj{q)t \ 
\ — u;{q) sin CO {q)t cosu;{q)t I 



- y e(--(^)-^/^)* f . I , ' ) , (5.2) 



so. 



s=±l 



where 



q') = g') + isu;{q)-'j{q, q'). (5.3) 

Thus, (4.9) becomes 

Kp{q, q') = l{27r)'' X\N^p{q, q') + N^^iq', q)) (5.4) 



with 



N{q, g') = ^ E / (e + A n ^^.(?^' • ( ? 

s •' \i=l ) i=\ ^ 

■^S3a;(g3)"^5(g3 + ^3)^54(^4, ^4) + ^^^{q^WsA^^^ q'sPiq^y^siq^ + 54)] (5.5) 





is4U!{q') 



where 



dfi = s(^q-^q^6 q'^ S{q' - q^) H dqidq'v (5-6) 

and wc replaced 2e by e. Wc got in (5.4) a factor {27rY for each lattice sum in the defintion 
of Mai3{q,q'), and a factor ^ for each of the four sums over Sj. Note that, because of the 
sum over s in (5.5), only terms that are even in s contribute, which means that the factor 

(j2i=i SiUj{qi) + itj gives rise to a delta function if the integrand is even in s, and a principal 



value if it is odd. 



We will look for solutions to (4.18)-(4.20) which are translation invariant in the directions 
orthogonal to the 1-direction. Thus, we look for solutions of the form 

G(g,g') = (27r)'^-^5(q + q')^(q,q',q) 

and 

G{x,y) = J e^(^^+'i'y)+^'i(^-y)^(q,q',q)(iq(iq'(iq. 
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where we write x = (x, x), and similarly for y. 

It will be convenient to change coordinates in the 1-direction. Let be a 27r periodic function 

on §Z2N X j^Z2N ■ Define 



on the set 



Then, 



^(P,k) = c/(p + k,p-k) 



^ (^)' ^^''^^'''^(q' = E (^) e^P(^+^)+*(^-^)^(p, k), (5.7) 



Indeed, each pair (q, q') gets counted exactly twice in the (p, k)-sum (because the pair (p, k) 
gives the same conribution as the pair (p + 2A'", k + 2A'"), with addition modulo 4A^), which 
accounts for the factor -\/2. Note that, in the N ^ oo limit, both sides tend to /[_7r,,r]2 ■(27r)"^ 
since the lattice spacing in the RHS of (5.7) is -\/2 • ■ 

Let p = (p, 0). We have then 

G{x,y) = I e'P^''+y^+'''^''-y^G{p,k)dpdk (5.8) 

where G{p, k) ~ g{p + k, p — k, k) and the integrals over the first components are Riemann 

sums, with the same convention as in (5.7). The function G is 27r periodic in all the variables 
and, moreover, 

G(p + 7f , A; - 7f) = G(p, k) (5.9) 

where it — (tt, 0). 

We will, from now on, work in the p, k variables, drop, in general, the tilde in G{p, A ), and 
we shall not distinguish between p and p, unless we need to stress that p is a number. Since 
the components of p (unlike those of k) other than the first one are always 0, this abuse of 
notation is harmless. 

Equations (5.5) then becomes 

S3C^(P3 + fcs) [a;(p3 + k^Y^6{2p^)WsSpA, k^) - cu(p4 + ^4)"^^ (2^4)1^.3(^3, k^)] (5.10) 

where 

du = 6{2p - Y^{pi + ki))6{^{pi - ki))5{p -k-p^- k4)dpdk, (5.11) 
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and k = {ki)f^i, p = {pi)i=i and where we used the identity (see (5.6)) : 

3 4 4 



q-J2Qi=P + k- {J2(Pi + h)-qi)=P + k- (^{Pi + ki) - q') 

1=1 i=l 1=1 

4 4 

= p + k- Y^ipi + ki)+p-k = 2p- Y,{Pi + ki) 



i=i i=i 

We can then write 

Afi2{p, -k) = ^(27r)3'^AV(p, k) , M2(p, k) = ^{27rf''X\n2ip, k) + n^ip, -k)) (5.12) 

with 

n{W)ip,k) E jnWsMM-P^)sMk^) ( , ) 

• 'uj{k^)-^5{2p^)WsM. kA - Pa) - (3 ^ 4)] i^,pk{dp dk) (5.13) 

where 

Uspkidp dk) = (^Siu{ki) + iej 6(2{p - Y^pi)^6{2p - ^ki)5{p - k - k^jdpdk (5.14) 

i i i 

In (5.13) we have shifted the fcj-integrals by —pi compared to (5.10). 
Let us intruduce the convenient notation 

uip, k) = {\{u{p + kf + u{p - kf))'^ (5.15) 

and 

5uj'^{p,k) ^uj{p + kf -uj{p-kf (5.16) 
Then the equations (4.18)-(4.20) read as follows in the p, k variables: 

u{p, k)^Q + i(( jr - r J)~ - M2(p, A;) - Ui2{p, -k)) = P (5.17) 

5u\p, k)Q + ( jr + r J)~ + M2(p, -A;) - J\fi2{p, A;) = (5.18) 

A;) J + (PF + rP)~ - N22{p. k) - C{p, k) ^ (5.19) 

We look for a solution that satisfies in x-space G{x,y) — G{—x,—y),Q{x,y) — Q{y,x) and 
J{x:y) = -Jiy.x) i.e. 

Qip, k) = g(p, -A;) = g(-p, A;) (5.20) 

J(p, k) = - J(p, -A;) = - J(-p, A;) (5.21) 



This is consistent since AfapiP-, k, W) = Mapi^—p-, —k, W{—-, — • 

We will see below that the lime^o^ is well defined for bounded W{p, k). Before that, let 
us outline our arguments. 
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6 Heuristics 



Before getting into the details we give an outline of the argument. Since the nonlinear term 
M depends only on Q and J, P can be solved from (5.17) in terms of them. Then (5.18) and 
(5.19) are two equations for the two unknown functions Q and J which we write as 

J)(g)+A^(Q,J)+A/'r(g,J)=(^) (6.1) 



where 

m,j)-{''^-'^"':^jJ'^f''''^) (6.2) 

and 

Nt{Q,J)^[ PrXli) (6-3) 

where P is expressed in terms of J and Q in (5.17). Nt collects the friction terms that have 
both linear and nonlinear contributions but will be treated as a perturbation. 

For large the sohition to this equation will be locally in the Xi coordinate a perturbation 
of the equilibrium state (4.22) given in the new variables by 

Qt{p, k) = Tu{p + k)-H{2^) (6.4) 

where 



5(p) = 27V5p,o. (6.5) 

Hence the importance to study the hnearization of M around that state. It turns out (see 
Section 9.2) that it is given by an operator which is a multiplier in the variable p: 

DN{Qt, 0) ( ^ ) (p, A;) = ( •] ) {k). (6.6) 
Cp is a matrix of operators 

£p=f^a ^:±n. (6.7) 



The translation invariant equilibrium (6.4) has support at p = (and due to the periodic- 
ity (5.20) at p = vr) . The nonequilibrium solution will also have most of its mass in the 
neighborhood of these points. Thus it is important to understand Cq. It will turn out that 

A,(0) = A,(7r) = 0, i^j (6.8) 

whereas £ii(0) is invertible. Invertibility of Cq would then follow from invertibility of £22(0). 
This, however, is not the case: £22(0) has two zero modes. 
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One of them is easy to understand. Eq. (6.4) is a one-parameter family of solutions to the 
7 = equations. Hence the derivative with respect to the parameter T is a zero eigenvalue 
eigenvector of the linearization £22(0) , i.e. 

£22(0)0;-' = 0. (6.9) 

There is, however, a second zero mode for £22(0) : 

£22(0)0;-' = 0. (6.10) 

While the first zero mode has to persist for the full Hopf equations due to the one parameter 
family of Gibbs states that solve them for 7 = 0, the second one is an artifact of the closure 
approximation. The nonlinear terms in the closure equations can be interpreted as describing 
phonon scattering and in our case only processes where two phonons scatter occur. These 
processes conserve phonon energy leading the the first zero mode and also phonon number, 
leading to the second one. The connected six-point correlation function which was neglected in 
the closure approximation would produce terms that violate phonon conservation and remove 
the second zero mode. However, for weak anharmonicity its eigenvalue would be close to zero 
and should be treated with some perturbation of the present analysis. 

The second zero mode leads one to expect that our equations have in the 7 = limit a 
two parameter family of stationary solutions which indeed is the case. These are given by (see 
Section 7.3) 

QtAp, k) = T(u;(p + kf - Au{p + k))-H{2p). (6.11) 

with J = and and P given by (4.18). For (6.11) to be a well defined covariance we need 
positivity of the denominator which holds ii A < m^. The zero-mode (6.10) is proportional to 
the derivative of Qt,a with respect to A, a.t A = 0. 

These considerations lead to the following ansatz for the solution 

Q{x,y) = Qo{x-y)+r{x,y), (6.12) 

where 

Qo{x - y) = Qt{x),a{x){x -y) + Qi{x,y) (6.13) 

The first term here is of local equilibrium form with slowly varying temperature and "chemical 
potential" profiles T(x) and A(x) (A, or more precisely T ■ A, can be considered as being 
related to a chemical potential, because it arises from the conservation of a number current, see 
section 7.2), and the second is a small perturbation (see (8.16) below and the remark following 
it). T{x) and A{x) are determined from the current conservation laws which arc projections 
of the equations (6.1) as follows. The operator £22(0) has two left zero modes. Projection 
onto them the equation (6.1) yields two nonlinear elliptic equations for the functions T'(x) and 
A{x) coupled to the rest of the variables i.e. J and r. For the latter taken in a subspace 
complementary to the zero modes the linear operator 600'^ + £p is invertible and J and r can 
be determined by fixed point arguments in a suitable Banach space as functionals of T(x) and 
A(x). The projected equations then allow us to determine the latter. 

Our main result, stated more precisely in Section 8, says that Q{x,y) is as above, while 
the currents corresponding to the two conservation laws discussed above and defined in section 
7, are linearly related, to leading order in |Ti — T2I, to the gradients of T(x) and A{x). T(x) 
is, to leading order in |Ti — T2I, linear in x, and therefore the currents are, to leading order. 
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7 Current conservation 



In the Introduction we recalled that the Hamiltonian structure of the dynamics leads to a local 
conservation law. We will show that our closure too has such a conservation law. 



7.1 Heat-current 

We start with a simple identity 

Lemma 7.1 / M22{p, k)dk = for all (bounded) W. 

Proof From (5.11) we see that / dkdvp^k is symmetric under the simultaneous interchanges 
S3 ^ S4, (ps, /cs) ^ {pi, ^4). Hence, the integral over k of the second term in (5.12) vanishes, 
because the integrand in (5.13) is antisymmetric and the rest of the integrand is symmetric 
under those simultaneous interchanges. □ 

Note that, by (5.8), 

Af22{x,x) = [ e^'P''Af22{P:k)dpdk (7.1) 



i.e., by Lemma 7.1, N'22{x,x) = 0, \/x. This follows also by inspection from (4.9) since 

{R{t)Go)^^ = -R{t)i2 (by (5.2)) and R{t)i2{z - y) = R{t),2{x -z)ifx = y. 

Consider now equation (4.20) restricted to the diagonal x = y. Let J' = ujJ + Juj. Then, 
uj'^J — JuP' = LoJ' — J'uj. Defining, for the unit vector in the /x-direction, 

= J'{x - ei„x), (7.2) 

we have, for cu given in (2.5): 

{u'^J - Ju''){x,x) = {-AJ' + J'A){x,x) 

= ~ ^M'^) ~ + ^M'^) + J'{x,x- e^) + J'{x,x + e^)) 

= 2 ^{J'{x, x + e^)- J'{x - e^, x)) = 2V • j (x), 

where (V^/)(x) = f[x + e^) — f{x). Note that, for other functions there will be also a 
current, but its form will depend on lo. In view of Lemma 7.1, (4.20), for x = y, reads 

V • J (x) + 7P(x, x){d^o + 5xjv) = 7(^1^x0 + ^a^xjv), 

or, since P{x, x) depends only on x, 

V • J (x) = 7((ri - P(0, 0))5xo + iT2 - PiN, N))5^n), (7.3) 

which is a current conservation equation; the heat current j has sources on the boundary. Since 
J depends only on x, we define j(x) = ji{x) and, so, 

V-j(x)=j(x+l)-j(x). (7.4) 
It is a useful exercise to rewrite this in the (p, A;)-variables. We have, see (2.5), 

uj(p + k) — uj(p — k) — 2 (^cos(p — k) — cos(p + k)^ = 4 sin p sin k 
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and so 



Su'^{p, k) — 4smpsmk(a;(p + k) + a;(p — k)). 



(7.5) 



Thus 




(7.6) 



where 




(7.7) 



for p G jf^2N- In (7.6) we used the fact that the k integral is 7r-periodic in p, due to (5.9), to 
write it as a function of 2p. It may be checked directly that this is the Fourier transform of 
ji{x) given by (7.2). 

7.2 Number current 

As stated in Section 6, the closure equations possess another, approximate, conservation law 
which we will derive now. Let p{p, k) be given by 



with u{p, k) given by (5.15), and project A22 now onto p instead of the function 1 as in Lemma 



Then 9 is 27r periodic. Unlike what happened in Lemma 7.1 ^ is not zero, but it will turn out 

to be very regular, sec Proposition 9.7 and Appendix B below, due to the fact that (7.8) at 
p = is a left zero eigenvector of the linearization of A/22- 

We will now integrate equation (5.19) multiplied by a linear combination of p and 1. For 
this, write it in the (p, k) representation. The covariance is 



p(p, k) = uj{p, k) 



(7.8) 



7.1. Let 



^(P) = / p{l,k)N22{l,k)dk=\{2^rX' j p[\,k)n,{lk)dk 



(7.9) 



C{p,k) = 2^(T, + T^e-^^''P) 



(7.10) 



and the friction term 



{VP + Pr)^(gi, 52) = 7 / dq [(1 + e*(^-''^)^)P(g, 92) + (1 + e^(''-''^)^)P(gi, q) 



So, after shifting q in the first integral by y and in the second by y, we get: 



{VP^PVy{p,k) ^jjdq [p(|,| + A;-p)+P(|,-|+p + A;)] (l + e^(^-2f)^). (7.11) 



Let 




(7.12) 
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with p given by (7.8). Integrating equation (5.19) multiplied by 77, we get 

/(x + 1) - /(x) - ^(x) + 7(x) = (7.13) 

where 

j'(p) = -i j dA;e-*P/2r7(p/2,A;)sink(cc^(p/2 + A;)+u;(p/2-A;))j(p/2,/c) (7.14) 

and 

7(p) = 7 / dkdvi^iv, k, q)P{l k) (1 + e^('i-P)^) (7.15) 

where 

^{p, k, q) = riip, k + '--l)+ri{p,k-'- + l) (7.16) 

is a smooth function. Note that the covariance dropped out from (7.13) because C{p,k) is 
independent of k, and we subtract from p its average in the definition of 77. Also, the latter 
does not contribute to 9 due to Lemma 7.1. 

Equation (7.13) is again a conservation law: 7 is a boundary term and, as we will see, 
9 vanishes in the limit of translation invariance. We call the current / the particle number 
current. 

7.3 Generalized Gibbs states 

We finish this section by checking that the states (6.11), with J — 0, are indeed solutions to 
the 7 = equations. Indeed, now, 

Ws{p,k-p) = QT,A{k)S{2p), 

with Qt,a = T{u{kf - Au{k))-\ So, from (5.13), we get: 

nip,k) = 5{2p)Y, j t{QTAki)sMk^) ( -^^J^^^) ) [oj{k^)-^QTAk^)-^{k^)~^QTAh) 

•(^ Siuj{ki) + iey^5{2p -^k,)S{p -k- ki)dk (7.17) 
Til is obviously even in /c, thus (5.18) holds. As for (5.19), write the [-] in (7.17) as 

-AT-^QTAk3)QTAki) [o{k^)-^ - uik^y-] , 
and use {x + ie)~^ — V — 2-nib{x) to get: 

n2(p,A;) = -27r5{2p)AT-''n2ip,k), (7.18) 

f ^ 

n2{p,k) = X) / n^^.^(^»)'S3S4 [c^(A:4) - coiks)] 

s 1 

4 

■d(j2siuj{ki))d{2p -Y.ki)5{p -k- k^)dk (7.19) 
1 
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(By s — > — s symmetry only the delta function contributes). The integral in (7.17) is supported 
on 

^SiLu{ki) = 0. (7.20) 
We will choose in (2.4) large enough so that (7.20) forces 

E^^ = (7-21) 

By symmetry, we may replace S3 by | J2f=i Si and, by (7.21) also by — |s4. Again, by symme- 
try, ssculks) may be replaced by ^J2f=iSiLj{ki) and by (7.20) by — |s4a;(/c4). Doing the first 
replacement for S3C<j(A;4) in the first term in the [•] in (7.19), and the second for s^uj{k^) in the 
second term, we see that (7.19) vanishes. Note that ^2 vanishes for all p and not only for 
2p = mod 27r. We shall need this later (see the derivation of (9.19) below). Hence (6.11) 
solves (5.17)-(5.19) if 7 = 0. 

Note however that these are not equihbrium Ti = T2 solutions, unless A — 0. This is 
because for A 7^ 0, the noise and friction terms in (4.20) will not any more balance each other. 

8 The space of local equilibrium solutions 

We will now describe the space where (4.18)-(4.20) are solved. To motivate our choice, consider 
the current conservation equation (7.3). Summing over x, we get, since j(x) is periodic: 



and so 



which is solved by 



Ti - P(0, 0) = -(T2 - P{N, N)) = jo (8.1) 



j(x + 1) - j(x) = jo{Sxo - So^n) (8.2) 



j{0)+jo xe[l,7V] 
J(x) = <! (8.3) 
j(0) xe[-7V + l,0] 



Since j(— x) = — j(x) we get j(0) = — |jo- The Fourier transform of this is 

f 1- 

.70—- 

i(p) = 



1 - e^P^ Ti ^ 

e»p — 1 W 



(8.4) 



p = 



(note that e*P^ takes only the values ±1). The current jo will turn out to be 0{1/N). 

We now describe the space to which functions such as J(p, k) belong. This space has to 
encode the 1/p singularity at origin in the equation (8.4) as well as the factor e*P^ coming from 
the fact that in x space there are two special points in the first coordinate, the origin and N . 
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Let Ti be the space of continuous functions /(p, k) on 

Q = {(p, A;) I p, k e ^'^an, P + k e ^Zsiv, k e [-tt, tt]'^-^} (8.5) 

that are 27r-periodic in all the variables, and invariant under (p, k) (p + 7r,k + 7r). We denote 
by ll/lloo the sup norm in Ti. 
Let 

n+ = {(p, A;) e Q I p e -^Zaiv} 

and r2_ = \ Then 
with 

/(P, ^) = ^)^+(2p) + /-(p, ^)a_(2p) (8.6) 

with 

a±(p) = l±e^P^ (8.7) 

and /± = i/|f7±. 

Define the discrete C"-norm for / G 'H± as 

||/||a= sup (|/(p,fc)| + |A|-"|/(p + Aei,A;)-/(p,^)| + H-'^|/(p,^ + //ei)-/(p,A;)|) (8.8) 
where X,Ijl E j^1-2N \ (note that p + A e Q± if p e Q± ). Set then, for f E H, 

\\f\\a^\\U\\a+\\f-\\a. 

Let 

r (e^P-1)-^ p^O 
d{v)-' = (8.9) 
i p = 0, 

and let, with some abuse of notation, 

rf(p)=e'P-L (8.10) 
We will then consider a space S of functions J(p, A;) of the form 

J = N-^5{2p)Jo + (A^(d(2p))-Vi + Ar"/2-i J2 + (Ard(2p))-3/V3 (8.11) 

where Ji(p, k) are in 7i. Define 

||J||5 = max{||J,|U,||J3||oc, i = 0,l,2}. (8.12) 

More properly, S is the space of 4-tuples ( Jq, ...,73) = J which is a Banach space with 
this norm. We identify J and J when no confusion arises. 
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Remarks. 1. Due to the friction and the noise terms (see (7.10) and (7.11)) we end up with 
the factors e^*P^. This leads to a factor e*P^ in functions such as j(p) and j'(p) . A Holder 
continuous function /, of exponent a, decays in x-space at least as Ixj""" whereas e*P^/ produces 
decay away from x = N i.e. |x — iV|~". Thus ||/||a < ^(1) gives rise to functions of x that are 
localized near the boundaries, x = or A^. 

2. The Fourier transform of (8.11) is as follows. The first term is constant in the 1-direction and 
0{1/N). The second one is also olO{l/N) because it involves a (discrete) Hilbert transform of 
a Holder continuous function. However, e.g. for the current it can produce opposite constant 
values for positive and negative x. The third term is also of this order whereas the second is 
e)(7V-i+"/2(| X |-" + I AT - X |-°)) i.e. C»(iV-i-"/2) for x far away from the boundaries. The 
J2 and J3 terms will be subleading corrections to the current and the profile. Note that no 
smoothness in the transverse variable k is assumed so that correlations will have a very slow 
decay in x space. This is due to the poor regularity properties of the "collision kernel" that 
enters the nonlinear terms, see Proposition 9.3. 

3. The J correlation function is odd in p and hence Jo = 0. However, with Q, we will encounter 
functions with a nonzero first term in (8.11). 

Finally, wc will extend these definitions to the functions like j(p) defined on jj~L2Ni the only 
difference being the occurence of o"-[;(2p) in (8.6) and (i(2p) in (8.11), instead of o?(p) here. We 
denote the space now by S. Hence j e 5' is of the form 



with Jo a constant. Of course, here, since j is defined on jj1-2n-, the space is finite dimensional 
and all norms are equivalent, but (8.13) is a convenient way to record the dependence on N of 
various terms. 

We will look for solutions to (4.18)-(4.20) with J & S and Q as follows. 

Q will have a leading "local equilibrium" term which we now describe. Recall that our 
equations have the 2-parameter family (6.11) of solutions when the friction vanishes. The 
leading term in Q will be of this form, where the constants T and A will be replaced by 
p-dependent functions. More precisely let T(p), A(p) be functions on j^Z2n, of the form 



i(p) = N-'5{p)jo + (A^d(p))-Vi(p) + N-^'-\ 



(8.13) 



T(p) = To5(p) + {d{-{p))-H{p) 
A(p) = Ao5{p) + (d(-(p))-ia(p) 



(8.14) 
(8.15) 



where 5(p) is defined in (6.5) and t,a & S. Define 



00 



2-n 



(8.16) 



n=0 



where * is the convolution and a;(p, k) is defined in (5.15). 



Remcirk Let Qt,a be the generahzed Gibbs state in (6.11), at p = 0, 




-1 



(8.17) 
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Comparing (8.16) with 



T{u^ - Auu)-^ = J2 TA'uj-^-'' (8.18) 

n=0 



it is not hard to show that the Fourier transform of (8.16) is 

Qo{x,y) = Qt{x),a(x){x -y) + Qi{x,y) (8.19) 
where (the Fourier transform of) Qi is in S, and, as it will turn out, 

Q^{x,y) = 0{T/N). (8.20) 
The first term is of a local (generalized) equilibrium form. 
We make now the following assumptions. We take 

^ = iV-iW4 (821) 

and large enough that 

^ SiU}{ki) = =^ ^ Si = 
(recall (7.20) and (7.21)). We shall use below the set of p's that are close to singularities: 

Eo = [-Po, Po] U [tt - po, TT + po] (8.22) 

with Po = BX^, where S is a number that will be chosen large below (see (10.16)). Then, we 
have the 

Theorem There exist Aq such that for < A < Aq, | Ti — T2I < r , with r = t(A), and for 
N > N{X) the equations (4.18)-(4.20) have a unique solution with J & S and 

Q^Qo + r 

with r eS. Qo is given by (8.14)-(8.16) with Tq = h{Ti+ T2) + C(AV), 

^(p)=^o(p)+^(p), (8.23) 

where, in x space, 

\{T, + T^) + t,{x) = + M(r2 - Ti), (8.24) 

and i & S, with 

\\i\\s^O{X'T). 

Moreover, Aq = C(t^), a & S, and 

\\a\\s^O{T'). 
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Finally, the currents j and j' 




where s(p) = d{—p)T * A{p), the conductivity matrix k is Holder continuous in p of exponent 
a, for some o; > 0, and is invertible and 0{X~^) for p & Eq. . 

Remarks 1. The shape of the profile, to leading order, is actually given by 



This can be seen, from the proofs, as follows: The full relation between J and t{p), which, in x 
space, is equal to — VT(,t), is given by (11.34) with Dp replaced by the full DN{qo) = C + C, 
see (9.22). The leading term in C comes from the n = term in (8.4) inserted in (9.21, 9.22). 
Going back to x-space, we get, within those approximations: T(a;)^j(a;) = — ^VT(a;), which, 
since j{x) is constant, imphes V(T'~^)(x) constant, i.e. (8.25). Of course, for r small, we can 
expand (8.25) in r and obtain (8.24), to leading order in r. 

2. The choice of constants. There are four parameters in this model: m^, the parameter 
in u!{k), see (2.5), A the strength of the nonlinearity, r that bounds the temperature difference 
and A^, the size of the system in the direction where there is a temperature difference. We 
choose enough so that (7.20) implies (7.21). Next we choose r small compared to A^, say 
0{X^) and we choose large so that we can use bounds like A^~" < r, for any a > 0. 

In the proofs, C or c will denote constants that can change from place to place. Wc also use 
in various proofs auxiliary functions denoted f,g, whose meaning is given in the proofs where 
they are used. 

Since A^ enters as a multiplication factor in the nonlinear terms, sec (5.12), it will be 
convenient to discuss in the next section the A independent nonlincaritics n(W), and introduce 
explicitely A^ in section 10 (because only part of the linear operator there is multiplied by A^) 
and in section 11, where A^ is used to make some Lipschitz constants less than 1. 

The value of a is determined by the degree of Holder continuity that one obtains in Propo- 
sition 9.3 below. We do not try to optimize that value (any a > suffices); in fact, it will 
convenient sometimes to assume that a is not too close to one, so that, e.g. the power on N in 
(11.21) is negative, and we shall imphcitely assume that. 

9 Nonlinear terms 

In this Section we study the nonlinear terms of our equations given in (5.12) and (5.13). Our 
goal is to show that their linearization defines a bounded operator on S and that the remaining 
nonlinearities define suitable Lipschitz functions on S, 

It is convenient to introduce a space for the functions T, A , analogous to S but stronger 
singularity at p = 0. We let E denote the pairs (To,t) with Tq G IR and t E S with to = in 
(8.13). They parametrize functions 



p{x) = T{x) 



-1 



(8.25) 



T{p)=To6{p) + d{-p)-H{p), 
(from now on, we shall identify p and p). We use the norm 

\\T\\e — \To\ + \\t\\s, 



(9.1) 



(9.2) 
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which has the convenient property, used often below, that ||<iT||5 < ||T||e. In this section, we 
shall work with (T, A) G Bg, where 

B, = {(T, A)eExE\ \\T\\e < C, \\A\\e < e} (9.3) 

where C is arbitrary and e is chosen small enough so that various series below converge. 
We expand (5.13) around Qq defined by (8.16). Let W — Qq + w 

n{W) - n{Qo) + Dn{Qo)w + h{w) (9.4) 

We discuss the three terms on the right hand side in turn. 

9.1 n{Qo) 

Consier first n((5o)- Let us start with a lemma, whose proof is given in Appendix B, and which 
controls convolutions between elements oi E, S and among themselves: 

Lemma 9.1. (a) Let T e E and j e S. Then T*jeS and 

||T*j||5<C||r|U||j||5 (9.5) 

(b) Let T,AeE. Then T * A e E and 

\\T*A\\e<C\\T\\e\\A\\e (9.6) 

(c) Let j, k & S. Then j' = j * k & S with 

\fo\ < CN-'\\j\\s\\k\\s, (9.7) 

UWa < CN-^j\\s\\k\\s, (9.8) 

IIj^IU < CN-'+-/'\\j\\s\\k\\s, (9.9) 

\j's{p)\ < CN-Hog{N\p\)\\j\\s\\k\\s, (9.10) 

and 

\\f\\s<CN-'+^/'\\j\\s\\k\\s (9.11) 

Using this Lemma, we get: 
Proposition 9.2 The function n(Qo) satisfies the following bounds, for (T, A) E B^. 

oo 

niiQo) = J2T*'* A*-{2p)g^{p, k)+m (9.12) 

n=l 

with Qn smooth functions hounded together with their derivatives (to any given order) by C" 
and m & S. Moreover 

\\m\\s<C{\\t\\s + \\A\\E) (9.13) 

\\ni{Qo){p,k)-n,{Qo){-p,k)\\s < C{\\t\\s + \\A\\e) (9.14) 

\\n2{Qo)\\s<C{\\t\\s+\\A\\E) (9.15) 
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Furthermore the functions on the LHS of (9.13), (9.14) o^i^d (9.15) are uniformly Lipschitz in 
T and A, for {T,A) e B,. 

Proof. From (5.13) and (8.16), we get 

2 

n{Qo){p,k) = / llT*A*''^{'2p^MP^,k^-p^)-'-''^ S^h) ( ) 

• ■a;(A;3)-'5(2p3)T*A*"^(2p4)c^(p4, h - P4)~'~"' - (3 ^ 4)] u,pk{dp dk) (9.16) 
with Uspk given by (5.14), n = (ni)^=i and uj{p, k) given in (5.15). We have, see (5.15), 

uj{p, k-p)= uj{k) + {e^'^' - l)C(l), (9.17) 

where the 0{p) term, is written in an unusual form, which records the fact that it vanishes also 
at p = TT, and which will be convenient later. Inserting this into (9.16), the leading term is 



EE / ^T*A*-^{2p,)uj{h)-^-''^s:,u{h 

n s ^ i=l 

■ [a;(A;3)-'5(2p3)T*A™^(2p4)a;(A;4)-'-"^ - 
We may do the pj-integrals to get 

oo .2 

= E^*'*^*"(2P)E E I[^{kr'-''^ssu;{k 

n=0 s T.ni=n i=l 

■ - (3 « 4)] v^t(dk) 

.|T-,x-(2rt(«;;fe^)) 



iS4!jj{k^ 

[3 ^ 4] Vspk{dp dk) 



is^^uj{ki) 



(9.18) 



where u^pk is like Vgpk above, but without 5{2(p — Y^Pi)), and where, in the last equality, we 
used the fact that the n = term has 71,3 = and therefore the [■] factor vanishes. Using (8.18), 
we see that fn{p, k) are the Taylor coefficients of the expansion in A of n2(p, /c), given by (7.18, 
7.19). Since 77,2 vanishes identically, we get 



fn{p.k)=Q. 



(9.19) 



For Qn we need to study rii in (7.17) in more detail. Proceeding as with ri2, we see that Qn 
are the Taylor coefficients of the expansion in powers of the constant A of 

n,{p,k) ^Y.j ilQo{ki)V {^Siu:{h)^ 

• 53(1 - uj{k^)uj{ki)-^\5 (2p - E^O ^(P-^- ^^)dk. 

Here, again because of the s — — s symmetry, only the principal value contributes. Consider 
first those terms with J2si — 0- We may replace S3 by — |s4 and S3a;(A;3) by 

1 3 1 1 ^ 

-Y,SiUj{ki) = --s^uj{k4) + -Y,SiUj{ki). 
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Thus these terms give 

Esi=0 1 

which is smooth in p. For the terms with Z^Si 7^ 0, V{-) has no singularity, (see (7.20), 7.21)), 
and they are smooth. Thus, the functions Qn are smooth and bounded together with their 
derivatives (to any given order) by C". 

The contribution to (9.16) of the {e'^'P - l)O(l) term in (9.17) is of form 

^ = EE / l[T*A*^^i2p,){e''P^ - l)r*A-3(2p3)5(2p4) 

n s i=l 

• G^{p,k,p,k)S(2p-2Y,Pi)^spk{dpdk) (9.20) 

where Gg is smooth. By an easy extension of Lemma 9.1 to convolutions of functions in E and 
S, each summand is in S and (9.13) follows easily. The series over n converge for H^H^; < e 
small enough. 

To get (9.14) we use the smoothness and periodicity to get 

gn(p,k)-gn(-p,k)^(e''P-l)0{l), 

and hence each summand is now in S. The bound (9.15) is obtained in the same way as the 
bound on m (by (9.19), the only contribution to n2{Qo) is similar to m ). □ 

9.2 Linearization 

Let us turn to the second term in (9.4) i.e. the linearization of n at Qo- We wish to separate a 
leading term when T2 — Ti is small. Therefore, we write, using (8.14) and (8.16), 

Qo = T+5{2p)uj-''{p, k) + Qo (9.21) 

with cu given by (5.15), and = |(Ti + T2) being the average temperature . Insert this into 
(5.13). The [ - ] term vanishes for the first term in (9.21), so 

Dn{Qo)w ^jCw + jC'w, (9.22) 

with 

iCw){p,k) = 2T,'Y. J {u;{k^)u;{k2)y^^u;{ks) ( ^^J^j^^^ ) 

[ioiksy^Ws^ip, k4 -p) - uj{ki)~'^Ws^{p, k^ - p) 

SiUj{ki) + ie) ~^ 5(2p-Y, h) S{p-k- k^jdk (9.23) 

(where we used a;(0, k) — uj{k), and i/;(p+7r, ki—p—ir) — w(p, ki—p), 1 — 3,4). Since k4—p — —k 
and k^ — p = k — ki — k2 (coming from ki + k2 + k^ ~ 2p — k^ — 2p — {p — k) — p + k), we see 
that £ is a multiplication operator in the p- variable 

{Cw){p,k)^{Cpw{Pr)){k) (9.24) 
28 



and Cp in turn is a sum of an operator which acts as a multiphcation operator in the subspace 
of even or of odd functions, and an integral operator: 

4 = + kp, (9.25) 

where, after shiflting kz by p, 

ik,v,)ik) = -2r„^.(p-^-'i: / '^3) jffy4l ( isjp-k) ) if'-*' 

with 

dfispk = [yI^ SiUj{ki) + S3'^(A;3 + p) + SiUj{p - k) + ie^ ^{^~^ j 11 (9.27) 

and (remember that — p = —k so that the operator acts as a multiphcation operator only 
on functions of definite parity): 

Mp{k) = -u;{p-k)\KpU;-\p + -)){k), (9.28) 

when it acts on even functions, and 

Mp{k)^u;{p-kf(KpU;-\p + -)){k), (9.29) 

when it acts on odd functions. Here, we integrated over k^ = p — k, and we used: 

4 3 

2p -^ki = 2p - "^ki - p + k, 
i=i i=i 

which, after shifting k^ by p equals k — ki. We shall discuss C in the next susection, but 
let now analyze further C 

Separating the real and imaginary parts of Wg (see (5.3)), one can also view the operator 
£p as a 2 X 2 matrix of operators the (J, Q) variables, where here Q denotes an arbitrary even 
function, and J an arbitrary odd one. We write it as: 

r ( J \- ( Ai(p)J + A2(p)g \ . Qm 
)-\C,,{p)J + tMQ ^^•''"^ 

this defines the operators Cij{p), i.j = 1,2. From (9.25), we see that each J0,ij{p) is a sum of a 
multiplication operator Mij{p) and a integral operator Kij{p). 

The operators here are acting on functions defined on Q, see (8.5). However, it is conve- 
nient to consider them as acting on functions defined on [— tt, tt]'^, which is always possible, by 
extending, say in a piecewise linear way, function on Q to functions on [— tt, tt]'^. With that 
identification, we may consider these operators to be acting, for all N, on the same space. 
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Moreover, a discrete Holder continuous function (see (8.8)) becomes, with such an extension, 
an ordinary Holder continuous function. The main property of these operators is: 

Proposition 9.3. Mp(k) is in p,k, for some a > 0, and the operators Kp are compact 
operators mapping C"{Qp{k)), where ^p{k) = {k\{p,k) G fl}, into itself. Moreover, in the 
norm of bounded operators on C°', they are uniformly bounded in N and in p, uniformly in 
N . Mp{k) converges as N —>■ oo to a function, while Kp converges to a compact operator 
mapping C"([—7r,7r]'^^ into itself. 

Obviously, this proposition implies that the operators Cij{p) define also bounded operators 
from S into itself. 

Proof. Looking at (9.26-9.27, 9.30), we see that Kp is a two by two matrix of integral operators 
whose kernel is a sum of terms of the form 

Ap{k,k') = A(j2sMki) + S30u{k' +p) + S40u{p-k)) 
■' i=i 

■5{k — ki — k2 — k')ps{ki, k2, k', k,p)dkidk2 (9.31) 

where A(x) = 5(x) or V Ps is C°° in all its arguments and we write k' for ks. 

Let us consider first A{x) = d{x), integrate over k2, and choose Si = +1, S2 = — 1 (all other 
terms can be treated similarly) . We obtain the integral 

j 5(uj{ki) — uj{k — ki — k') + s^uj{k' +p) + SiUj{p — k)^ps{ki, k — ki — k', k', k,p)dki. (9.32) 

Now, by (5.1), / dki is actually a discrete sum over the first component of ki and an integral 
over the last two components. Fix a value of k}, replace k' — k hy k' , write ki — — | + g2, 
k\ = —\ + qz and use lower indices, k[,k2,k'^ for the components of k'. We get, using the 
explicit formula (2.5) for cu, 

(9.32) = j d(smq2 - sin(g2 - k'^) + sin 53 - sin(g3 - ^3) + f{k2, k'^, X))psiq2, Qs, X)dq2dq^9.33) 

where we write A for the set of variables {kl, k[, k,p), the integral is an ordinary, not discrete, 
one and 

7(^2, fcg. A) = sin qi — sin(gi — k[) + Szijj{k' + p) + SiU}{p — k) (9.34) 

We shall now study the singularities of (9.33) in /c2, k'^ and the smoothness of (9.33) in A 
away from the singularities. For notational simplicity, we set ps = 1; since ps is smooth, this 
does not affect our arguments. Let us denote by Q the argument of the delta function and shift 
the integration variables: let q^ — ri-i + where = ^k^. Then 

2 2 
Q = ^(sin(ri + yi) - sm{ri - yi)) + f {2y , X) = 2 ^ cosn sin + /(2y, A). (9.35) 

1=1 i=i 
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Next, change variables to 2 — 2cosri = x'^. Our integral becomes 



f 

I{y,X)= / S{xlsmyi + xlsmy2- g{y,X))Y[h{xi)dxi (9.36) 

with h{x) = (1 - x\ e [0,2] and 

g{y, A)) = /(2y, A) + 2(sinyi + sin y^). (9.37) 
For yi^ Q and 5^ 7^ / is bounded by 

|7(y, A)| < C| sinyi siny2r^(l + | log \g{y, A)||) (9.38) 

(the log term is absent when yiy2 > 0). From (9.34) and (9.37) we have 

g{y, A) = 2s3(cos(2yi) + cos(2y2) + 2(sin|/i + sin 7/2) + 5''(A) 

with g' smooth. Thus the bound (9.38) is intcgrablc in y, uniformly in A. Moreover, the 
Holder derivative of order a in A also remains integrable for a small enough. Hence for all 
bounded functions /, / Ap{k, k')f{k')dk' is in A = {kl, k[, k,p). This in turn implies that 
/ Ap{k, k')f{k')dk'dk\ (where now the integral includes a Riemann sum over k[, kl, is C°' in 
k,p, since each term in the Riemann sum is in k,p. This means that each matrix element of 
Kp maps bounded functions into Holder continuous ones. Moreover, all the bounds are uniform 
in A^, since the bound (9.38) is independent of N, and taking the Riemann sums preserves that 
property. 

To obtain compactness, let a' the Holder exponent obtained above. Using the Arzela- 
Ascoli's theorem, one easily shows that the unit ball in is compactly embedded in for 
a' > a. Since Ap is bounded from into itself and maps C" into C"' it a compact operator 
from into itself . 

Obviously / Ap{k, k')dk' is also in p, /c, so that each matrix element of Mp{k) is in 
p,k. 

Next, we get: 



\Apf{k)-Ap,f{k)-{Apf{k')-Ap,f{k'))\ < C||/|Umin(|A;-A;r',b 



p 



/|a'^ 



< C\\fUk-kr'^'\p-pT'^^ (9.39) 

so that, choosing a — a'/ 2, we get: 

\\Ap-Ap,\\^<C\p-p'\'' (9.40) 

where || ■ \\a is the operator norm on bounded operators of into itself. 

Finally, it is easy to see that the Riemann sum of a function converges to the corre- 
sponding integral, with, for the sum in (5.1), an error (!?(iV~"). This then implies the claims 
on convergence as A?" — > 00 made in the Proposition. □ 

In Appendix A, we will extend this result. 
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9.3 Nonlinearities 

Let us now turn to the rij's defined in (5.13). They are linear combinations of functionals of 
the following form 

3 3 

u{fi, /2, f3){p, k) = fell Mp„ h - Pi)5{2p - ^ 2p,)dMdk) (9.41) 

i=i 1 

2 3 

v{fi,f2,f3){p,k) = fGl[Mp,,k,-p,)h{ps,p-k-p,)6{2p-J2^p,)dMdk) (9.42) 

i=i 1 

where G is some smooth function of all the variables (defined in the continuum {N — oo) and 
restricted to the discrete) and 

3 -1 3 

IJ'P,k{dk) — f y^ SiUj{ki) + SiUj{p — k) -\- iej k^ — p — k)dkidk2dk^. (9.43) 

1 1 

Indeed, (9.41) corresponds to the second term in the bracket in (5.13), after integrating over 
/c4, P4, while (9.42) corresponds to the first term, after integrating over /C4 (i.e. replacing k^^ by 
p — k) and ps, and relabelling ^4 as p^. Since Q = Qq + r , r & S and Qq is given by the series 
(8.16) we only need to consider /j(p, k) — F{2p) ior F e E or /j e S. We have then the 

Proposition 9.4. Let m be of the form (9.4I) or (9.42). Then 

(a) Let fi{p, k) = Fi{2p) with Fi e E. Then 

||m(/i,/2,/3)||E<Cni|F,|U (9.44) 

i=l 

(b) Let fi e S, f){p, k) = Fj{2p) with Fj eEforj^i. Then 

IM/i,/2,/3)||5<C^||/.||5nil^ilU (9-45) 

(c) Let fi, fj eS, fte E. Then 

||m(/i,/2,/3)||5 < C N-'^+^fkhUllMs (9.46) 

l^k 

For the proof, see Appendix B. 

The operators C in (9.22) are not multiplication operators in p. Inserting the expansion 
(8.16) in (5.13), we see that Cw is a sum of terms of the form discussed in Proposition 9.4. (b) 
with Fi of form T*>1*" with n > or d{-2p)-H{2p) . Proposition 9.4.(b) then gives the 

Proposition 9.5. The operator C : S ^ S is bounded, for (T, A) e Be, in operator norm by 

m<Cms + \\A\\E), (9.47) 
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and it is uniformly Lipschitz in T and A for (T, A) E B^. 
Proposition 9.4. (c) gives immediately for the term n{w) in (9.4) : 

Proposition 9.6. For \\w\\s < 0{1), and {T,A) e B, \\n{w)\\s < CN-^+"\\w\\l. and it is 
Lipschitz in w, T and A with constant CN~2+°'. 

We still need to discuss the function 9 given by (7.9). Its main property, proven in Appendix 
B, is: 

Proposition 9.7. d'^O is in S, for {T,A) e B^, with, for \\w\\s < 0{1), 

\\d-'e\\s < CX'iWtWs + \\A\\e + \\r\\s\\J\\s + II Jill). (9.48) 

and it is uniformly Lipschitz in T and A for (T, A) G B^, and in J and r, with constants 
CX'^ .Moreover, 

\e{0)\<CN-\ (9.49) 

10 Solution of the linear problem 

In order to solve equations (5.18), (5.19) we need to study the invertibility of the linear operator 

5u;'{p,k)^[ 1)+jC, (10.1) 

where 600"^ {p, k) — uj^{p + k) — u}^{p — k) and Cp denotes the linearization around the first term 
in (9.21) of the nonlinear terms in (5.18, 5.19). Cp is given explicitely by adding or subtracting 
to (9.23) a term with k -k, see (5.18, 5.12) , and multiplying it by |(27r)^'^A^, see (5.12). 
The fact that C in (9.22) is a small perturbation of C follows from Proposition 9.5, for r small 
enough, since, as we shall see in the next section, we shall solve our equations in a space where 

the RHS of (9.47) is of order r. So, to invert Suj^{p,k) ^ ^ ) DN{Qq), it is enough to 

concentrate our attention to (10.1), where Cp is written as a two by two matrix as in (9.30). 

Because of the zero modes, in order for (10.1) to be invertible, it needs to be restricted to the 
orthogonal complement of the zero modes (which occur at p = 0, tt). Let T-Lp be the Hilbert space 

L'^{j;^l_2N X [—T^i ^(p, kydk^ and let be the projection to the orthogonal complement 

of |a;(p, A;)~^,a;(p, A;)~^| in Hp, and P — 1 — P^. Note that the scalar product of f{k) with 

u!{p,k)~^ equals / f{k)dk, which implies 

P^l = 0. (10.2) 
We shall study the operator Dp, defined as: 

Vp = U (duj^ai + Cp) n, (10.3) 
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where we use the shorthand notations "^i — ^ ^ q ^ : 



^0 y 

for p e Eo, with Eo defined in (8.22). 

Vp = [duj^ai + Cp) , (10.4) 

for p ^ Eo. 

We shall make the following assumptions, that we shall verify later for the operator Cp : 
1. A,(0)=A,(7r) = 0, i^j (10.5) 

2. 3 c > such that 

£ii(0) < -cA^ (10.6) 
P^C22{0)P^ > cX^ (10.7) 

where the inequalities holds for operators in (^-^^2N x [— tt, 7r]'^~^,uj'^ {k)dk) restricted to func- 
tions that are odd in k (for Cu) or even in k (for £22)- The same inequalities hold for £11(77), 

>C22(7r). 

3. 3 C < 00, independent of N, such that, \/p,p',\/i,j = 1, 2, 

\\C,j{p) - C,j{p')\\ < C\p - pT (10.8) 

where || ■ || is the operator norm of bounded operators mapping (^j^Z2n x [~7r, tt]'^"^) into 
itself. 

4. 3 c> such that, Vp, V/c, V/ e R^^ 

(5a;Vi + M(p,A;))/| > c(A^ + | sinp| \ smk\ )|/| (10.9) 

where I/I = |/i| + |/2| 

5. The kernels Kij{k, •) e L^^'^{j^Z2N x [-tt, tt]'^ ^), for some > 0, and a norm C(A^). 

Proposition 10.1 Under assumptions (1-5) above, for any B in (8.22), 3Ao such that, for 
A < Ao and for all p e -§^~1-2Ni T^p is invertible, and 3C < 00 such that VA < Aq, 



< (10.10) 



where the norm is the one of operators in C" {j^'Z2n x [— tt, tiY^^^ © C" (^;^Z2Ar x [— tt, tt]'^^^ 
Proof. We consider |p| < |. For |p| > | the proof below can be repeated with replaced by 

TT 
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Let us first consider |p| < i?A^, where tlie constant B, independent of A, will be specified 
later. Write p — aA^, with \a\ < B. We have, for \p\ < B\^, 

Su\k,p) = a\^^{k) + 0{X*) (10.11) 

where we expand Suj'^{k,p) (which vanishes at p = 0) in p, the first term in (10.11) is the one 
hnear in p, with (p an odd function of k alone, and the second is C(p^) for \p\ < B\^. Moreover, 
due to (10.8), 

A,(p) = £,,(0) + O(A2+") (10.12) 
where 0(A^+") is a bound on the operator norm {Cij{p) has a factor A^). So, we have 

Vi + £p = A^ (a^{k)ai + Co) + 0{\^+'') (10.13) 

with Cq — Cq/\^, which is A-independent. Thus if we show that 3 C{B) such that Va, \a\ < B, 

II \a(a^{k)ai+ Co) n]"^ || < C{B) < oo, (10.14) 

we obtain, from (10.13) and a resolvent expansion the bound (10.10) for |p| < BX'^, provided 
that A is small enough, given B. 

To prove (10.14), we observe that the spectrum of the multiplication part of aip{k) ai + Cq, 
i.e. aip{k)(Ji + MSMl outside a ball of fixed radius around zero. This follows from (10.9), 
using the approximations (10.11) and (10.12) for Soj"^ and M(p, k). Hence, since adding to this 
a compact operator docs not change the essential spectrum (note that the projection operator 
P-*- adds a rank 2 operator), (10.14) will hold, provided that we show that aui + £(0) has no 
zero eigenvalue. But solving for /2 the first of the two equations 

/ £n(0) am \ h{k) \ 

(we use the fact that, by (10.5), Aj(0) = 0,i 7^ j), substituting in the second equation, and 
taking a scalar product in L^^-^Z2jv x [— tt, tt]'^"^, a;^(A;)(iA;j, with /i, we get: 

(P^/i,£n(0)P^/i) - a^(P^/i,(^At(0)<^PVi) = 

which is impossible (for all a's) because of (10.6) and (10.7). This finishes the proof of (10.14) 
and therefore of (10.10) for |p| < PA^ 

Consider now |p| > PA^. Since adding a compact operator does not change the essential 
spectrum, which is bounded away from zero by (10.9), it is enough to show that the equation 

{5uj^ar + C,)f = iif (10.15) 

with / = ( /^(A;) )' '^^^ ^^^^ solutions for |//| < cA^, and c > 0. Let E — 

{/c||sin/c| < h] for some constant h to be specified later. Then, from (10.9), (10.15) and the 
fact that Cp — Mp + Kp, we get, for k ^ E: 

c'BbXy{k)\ < \{Su;'a, + Mp)f{k)\ < \\Kp\\ ||/||oo + |/^| 
< (CA' + 1/^1 
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using the fact (see the proof of Proposition 9.3) that Kp maps bounded functions into C" 
ones which, in particular, are bounded, and that \\Kp\\ < C\^. Given /3>0,c>0, c'>0 
C < oo and 6 > 0, we can choose B — B{I3, c, c', C, 6), independent of A, so that this imphes, 
if < cA^, 



1/(^)1 </3||/l|oo 

for k ^ E. 

Consider now k e E. Then, we get from (10.9), (10.15) : 

cX'\f(k)\<\(5L0^a, + Mp)f(k)\ 



(10.16) 



< , 
' /fc'eE 

We bound, using (10.16), 



K{k,k')f{k') + I K{k,k')f{k')dk' 

Jk'^E 



[ K(k,k')f(k')dk' 

Jk'^E 



<CX'f3\\fl 



(10.17) 



(10.18) 



and since A;' G E means, for b small, that k' must be close to zero or close to tt, we get, using 
assumption 5 and Holder's inequality 

/ K{k,k')f{k')dk' < \\K{k,-)\\i+,{cby>/'+'^>\\f\\o. < CAV/^+^II/IIoo. (10.19) 

Jk'eE 

Inserting (10.18) and (10.19) in (10.17), we get |/(A;)| < (C"/9 + ^) ||/||oo, for k e E,hj 
choosing b small enough. This, combined with (10.16) implies f — 0, for /3 small enough, if, 
say, < Thus, there is no non-zero solution of (10.15) in that ball and (10.10) holds. 

□ 

Let us now check that the operator Cp has the properties (1-5) above. 

To do that, we must first write explicitely the operators £.^(0), Cij{7i). Let us start with 
p = 0. Using (9.23, 9.30), we get, since only terms that are even in s give a non-zero contribution 
(see (5.5) and comments afterwards): 

jCu{0)Q{k) = 2T^Y1 I {co{ki)u;{k2)y's3u;{ks) [u;{k)-^Q{ks) - cuik^r^Qi-k) 

■v{{j2siUj{K)y'y{^k)5{k + k^)dk- {k ^ -k), (10.20) 

—k) term comes from (5.18); (10.20) vanishes because Q and u; are even in k, 



where the {k 
see (5.20). 



Uj{k3 



-ujiksY 



LuJik) = 2T^Y1 I (^{kiMk2)) ^s^uiks) 

s 

■6(j2sMki))s{J2ki)^ik + kA)dk-{k -k) 
which equals twice the first term, since J is odd in k (see (5.21)) and uj even. 



uj{k) 



;io.2i) 
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>C2i(0) J vanishes by symmetry, like (10.20) (there is a + (A; — > —k) term in (5.12) and J is 
odd). This proves property 1, for p = 0. 
To prove point 2, write: 



(10.22) 



{C22Q){k) = 4T^Y1 / [^{kiMk2)) ssuj{ks)s4Uj{k4) [uj{k4)-^Q{ks) - uj{k3)-^Q{k) 
■5{^SiUj{ki))5{Y^ ki)S{k + k4)dk. 

Taking the scalar product of (10.22) with Q{k) in L2(f Zsat x [-7r,n]'^-^,u;^{k)dk), and 

/I 3 \ 

replacing in the second term in [ — ], 530; (/cs) by I - ^^SiUj{ki) J , using symmetry, and then by 

i=i J 

-S4Uj{k4 ) , using the delta function, we get: 



s i=l 

■5(j2siu{ki))5(j2 ki)S{k + k4)dk. 



S3u;^{k3)s4Uj\k4)Q{ks)Q{k) + ^u;%k4)\Q{k) ' 



(10.23) 



Using k — —k4, the symmetry between the kiS and the evenness of cu, Q, we can write 
(10.23) as: 

?E / li^iki)"^ \Y,sMhfQ{ki)\'s(Y,sMki))s{Y.ki)S{k + k4)dk. (10.24) 
A similar computation, starting with (10.21), leads to 

(J,£nJ) = -'^Yl [ nc^(^0"'|E'^(^0'-^(^i)f <^(E«i^(^0)<^(E^O'^(^ + ^4)rfMlO-25) 

s i=l 

To conclude the proof of (10.6), (10.7), we need the following 
Lemma 10.2. Let f be a Holder continuous function from T*^ into R, with d> 3, satisfying 

f{h) + f{k2)=f{ks)+f{k4) 

on the set 

{{ki)t=i\ ki e T'^, i = 1,...,4, ki + k2 = k^ + ki, uj{ki)+u{k2) = uiks) + Lu{k4)}, 

then, 

f(k) = au;(k) + b. (10.26) 



Now, consider first (10.24) and (10.25) in the N 
replaced by an integral. 



00 limit, i.e. with the sum in (5.1) 
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To apply the Lemma, use the fact that, by (5.21) and relabelhng indices, we may assume 
that Si = S2 = 1, S3 — S4 — —1 and change ^3 — > —ks, — > —k^. Then, the Lemma apphed 
to / = uj{kYQ{k) or f — uj{kyj{k) imphes that (10.24) and (10.25) cannot equal zero unless 
J = 0, since J{k) is odd, or unless Q{k) = auj{k)^'^ + huj(k)^^ which, for Q = P^Q, i.e. 
for Q orthogonal to a;(/c)~^ and to u!{k)~^ implies a = b = 0. Since each Cn{0) is the sum 
of a multiplication operator which is bounded away from zero (see (10.28), (10.29) below for 
p = 0) and a compact operator, either there is a zero eigenvalue or inequalities (10.6)-(10.7) 
hold. Since we showed that there cannot be a zero eigenvalue (10.6)-(10.7) hold for p = 0, and 
= 00. 

To conclude the proof for large, but finite, we simply use the convergence results stated 
in Proposition 9.3, which imply the convergence of (10.24) and (10.25) to their iV = cx) limit. 

The statements for p = tt follow from (9.23) and the fact that ^(Tr, k — tt) = w{p, k). 

The bound (10.8) follows from Proposition 9.3, and point 5 follows easily from the bounds 
(9.38) on the kernels Kij{p){k, ■), given by (9.31) . 

Finally (10.9) holds because, for given p, k, (Suj'^ai + M{p, k)^ is a 2 x 2 matrix, and the 

lower bound (10.9) holds if the eigenvalues of that matrix satisfy 

|A*i(p, ^)| > c(A^ + I sinpll sin A;|), i = 1, 2. 

To prove this lower bound, it is enough to prove it for the square root of the absolute value 
of the determinant of the matrix 5cu^ai + M{p, k), which equals 

Mn{k,p)M22{k,p) - (5u^ + Mu{k,p)){Su;^ + M2i{k,p)) (10.27) 

We can check, from the explicit formulas (10.21), (10.22), in which the multiphcation operator 
corresponds to the last term in the [ — ], that, for all p, 

Mnik,p) < -cA^ (10.28) 
M22{k,p) > cX^ (10.29) 

for c > 0. The signs and the factor are obvious, and to get a non zero contribution, we need 

only to check that J2siU!{ki) vanishes for some ki, k2, k^, k4^, with the constraints k^ = p — k, 
3 

kj — p + k. Choosing si — +1, S2 — —2, S3 = +1, S4 = —1 (which can always be obtained 

i=l 

by relabelling indices), and inserting the constraints, this means 

u!{ki) — uj{k2) + u}{p + k — ki — k2) — a;(p — k) 

which vanishes for ki = k2 = k. 

Now, since Mij{k,p) for i ^ j vanishes at p = or tt and at A; = or tt, we get, using 
Proposition 9.3, 

\Mij(k,p)\ <CA2(|sinA;| |sinp|)"^^ (10.30) 

for i 7^ j. 

Inserting (10.28), (10.29), (10.30) into (10.27), using \Suj^\ = 4|sinA;| |sinp|, we get that 
^ijj ^ 7^ J is small compared to if |sin/c| | sinp| > c'A^, for c' small (in which case. 
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both terms in (10.27) are negative), and that [Suj'^ + Mij(k,p)j is small compared to A^, i.e. 
compared to Ma, otherwise. 
We are left with the 

Proof of Lemma 10.2. The proof follows closely the one of Proposition 12.1 in [29], which 
itself is inspired by [8] . 

Let us first assume that / is twice continuously different iable. The hypothesis of the Lemma 
imply that f{k) + f{k') is constant \/k, k' e T'^, with k + k' constant and u!{k) +u!{k') constant. 
Therefore, there exists : IR x T'' ^ R, of class C^, such that 

f{k) + fik')=g{u;ik)+uj{k'),k + k'). 

Writing k = {k")i=i, A;" e T, u; = u}{k'),u}' = u}{k'), we get 

daf{k) = d^g{uj + uj' ,k + k')daio + dag{uj + uj',k + k'), 
daf{k') = d^g{u} + u}',k + k'), daou' + dag{oJ -\-uj\k-\- k'), 

where da — gfj, = Subtracting these two equations, we get 

{dj{k) - dJik')) = d^g{u + u;',k + k') (OMk) - dMk')) (10.31) 

Multiplying first (10.31) by (j)f^uj{k) — dpuj{k')^^ then rewriting the resulting equation by 
exchanging a and j3 wc get, for all a, (3: 

{dafik) - dj{k'))(dpuj{k) - dpcuik')) = {df,f{k) - dpf{k'))(dauj{k) - dacuik')). 

If we differentiate this identity with respect to k^, we get: 

dad,fik){d0uj{k) - dpuik!)) + (dj{k) - dj{k'))dpd,uj{k) 
= dpd,f{k)(d^uj{k) - d^k')) + (dfffik) - dpf{k'))dad^u;{k). 
Differentiating now this with respect to A;^, wc get: 

d^d^f{k)dpdsLo{k') + d^dsf{k')df,d^u;{k) = dpd^f{k)d^dsuj{k') + dpdsf{k')d^d,.Lu{k). (10.32) 

Now, for uj{k) as in (2.5), we have dadf3Uj{k) = Sais cos k^. Using this and choosing in (10.32) 
a ^ 'J 7^ f'^ = S, we get dad^f{k) = for a 7^ 7. This holds first on a dense set ky 7^ ±| 
(cos k^ 7^ 0) and then, by continuity, everywhere on T''. Then, choosing a = 7 7^ /5 = (5, we get 

dlmdlu;{k') = d}f{k')dluj{k). 

which implies that d'^f{k) — ad^uiik) for a constant a. Integrating, we get 

f{k) = auj{k) + h + ck, 

for a,b,c & R, and we get c = from the fact that / is a continuous function on T*^. 

This finishes the proof of / of class C^. For / merely Holder continuous, we interpret all 
the (linear) identities above, and all the derivatives in the sense of distributions, and we obtain 
the same conclusion. □ 

Remark. The Lemma holds by assuming only that / is a distribution, but we do not need 
this. It is crucial here that the dimension d > 3 in order to be able to choose a 7^ 7 7^ For 
counterexamples in d = 1, see [19]. 
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11 Proof of the Theorem 



In this section we solve equation (6.1). Due to the presence of zero modes in the operator 
£22(0) we need to consider separately (6.1) in the complement of these zero modes, at least for 
p & Eq, where 

Eo = [-po,Po] U [vr -po,7r + po], 

with pq = BX^, was defined in (8.22), and the projection of (6.1) onto the zero modes. The 
solution to the complementary equations leads to the Fourier Law, i.e. an expression of the 
currents {j.j') in terms of the temperature T and chemical potential A. The solution to the 
projected equation determines finally T and A. 
We look for a solution of the form 

W ^Qo + w, 

where Qo = Qq{T, A), for some functions T, A, is given by (8.16), and w, also written as a pair 
(J, r), is as follows. Let 

Ws^w x{p e Eo) 

Let P be the projection in L'^(u}{p, k)'^dkj to the span of |a;(p, A;)~^, uj{p, A;)~^| and P-*- — 1 — P. 
We demand 

P^rsip,-) = Tsipr), P^Eq 

Given a function / that is Holder continuous in p Eq, let / denote a linear extension of 
/ to jf^iN- We have ||/||q, < C||/||q;. We proceed similarly with elements of S or of E. Now 
write 

w = Wg + We 

where Ws is the extension of Ws defined above. The function satisfies Wj^ip) = for p e £^o- 

Since from (7.10, 10.2) we have PC — C, eq. (6.1) can be written, for p e £'0, as a pair of 
equations: 



n 



6u;^ 
Su^ 



( Q)+J^{Q,J)+J^r{Q,J)j=o (n.i) 



where 



A) 

and 

P{5uj^J + Af2+Afr2) = PC. (11.2) 
For p ^ Eq, we will solve directly (6.1): 











j 


(^) 
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We look for solutions where {T,A), defined for p e Eq, is such that {T,A) e Bi, where 
BiGE X Eis the ball 

\t - T^5{p)\l < B,T (11.4) 

\\A\\e < B^r^ (11.5) 

with 

r+ = |(ri + r2) (ii.e) 

and Bi will be fixed later to be 0{1). For r small enough, we have Bi C Be, defined in (9.3), 
so that the estimates of section 9 can be used. As for w, we look for Wg with Ws in the ball 
B2C S ®S, given by 

^2-{(J,r) I UJ,r)-iJo,fo)\\s<B2T} (11.7) 

where (Jo)^o) is given in (11.36) and bounded in (11.37). Finally, we shall choose e B3, 
where 

Bs^{{J,r) \\\J\\s+\\r\\s<B^X-\ Jp ^ Vp ^ for p e Eo}. (11.8) 
Our Theorem is a consequence of the Proposition below and the Remark following it. 
Proposition 11.1. Let \, t, N be as in the Theorem. 

(a) Given {T,A) defined for p G Eq, such that {T,A) G Bi, and given W( G B3, there exists 
a unique Wg, such that Wg G B2 and such that Qo{T,A) + Wg + wg solves (11.1) for p G Eq, 
Moreover, Wg is Lipschitz in {T,A,Wi) with Lipschitz constant O(A^). 

(b) Given wi G B3, there exists a unique {T,A) defined for p G £^0? such that {T,A) G Bi, 
and such that Qo{T, A) + Wg + solves (11.2), for p G Eq, where Wg{T, A, we) is the solution 
obtained in (a). The pair {T,A) is Lipschitz in we with Lipschitz constant 0{t). 

(c) There exists a unique W£ G B3 such that Qq{T,A) + Wg(T, A,Wi) + solves (11.3), for 
p ^ Eq, where {T,A) — {T,A){wi) is the solution obtained in (b). 

Remcirk. Moreover, the precise bounds stated in the Theorem will be given in the course of 
the proof: see (11.81), (11.86), for the statements about t, (11-5) for the ones on A, and, for 
the currents, see (11.7) the Remark at the end of subsection 11.1, specially (11.43). 



11.1 Fourier's Law 

Let us prove first part (a) of Proposition 11.1. For simplicity of notation, we shall write Qq 
for Qo(T, A), and drop the tilde on T, A. Recall that Q — Qq + r. The leading inhomogeneous 
term in (11.1) is —Slu^Qq. Using (8.16), we write: 

Su^Qq = t(2p)pi(p, k) + d{-2p){T * A){2p)p2{p, k) + p3(p, k) (11.9) 
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where t{p) = d{—p)T{p), 



Pi{p,k) = d{-2p)-^Suj^{p,k)uj{p,k)-^ (11.10) 
P2{p,k) = d{-2p)-'du;''{p,k)u;{p,k)-^ (11.11) 



are smooth functions and 

oo 

ps{p, k) = 5uj\p, k)Y.T* A*Mp. ^)"'"" 

n=2 

is in S with 

IIP3||5 < C\\A\\l. (11.12) 

The nonhnear term J\f{Q, J) in (11.1, 11.3) was studied in Section 9. It is given by (see eq. 
(5.12, 5.18, 5.19)) 

N{p,k)--{27r) A + J (11.13) 

and n is the sum (9.4). From Proposition 9.2, we infer 

mQo,0)\\s<CX'ms + \\A\\E). 
From (9.22), (9.24), the definition of Cp in section 10 and Proposition 9.5, we get 

WDMiQo, 0) - Cp\\ < CX\\\t\\s + \\A\\e) (11.14) 
where on the LHS the norm is the operator norm in <S © <S. Finally, Proposition 9.6 gives 

\\J\f-J\f{Qo,Q) - DAf{Qo,Q)Urf\\s < CN--2+-{\\J\\s + \\r\\sf. (11.15) 
So, combining those estimates, we get that , for T, A, J, r, as in the theorem: 

\\M{p, k) - Cp{J,r f\\s < CW (11.16) 

and is Lipschitz as a function of T, A with constant CA^, and as a function of J, r, with a 
constant CN~2+°'. Consider next the function Ar, defined in (6.3): 

jVr = (F J + JF, FP + PF)^. (11-17) 

Let us start with FJ+ JF, given by (7.11) with P replaced by J. RecaUing the definition (8.11), 
using 

J \d{p)\-^dp < Clog TV, (11.18) 
that follows from |(i(p)|~^ < C|p|~^, ioi p ^ and \p\ > cN~^, we get 

J \Jiq,q + k-p)\dq<CN-'+-/^J\\s, (11.19) 
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uniformly in k,p, since the singularities in (8.11) affect only the first variable and log < iV"/^. 
Now, use the fact that, for functions on the 7r/2N lattice. 



<CN"\\f\U (11.20) 
and identify F J + JF as an element of S of the form F J + JF = (0, 0, 0), to get: 

\\TJ+JT\\s<C^N'^N'-'^/'\\ J \J{q,q + k-p)\dq\\^ 

< C-fN''\\J\\s < CiV-'+^"/^|| J||5 (11.21) 

where in the last inequality, we use the definition (8.21) of 7. F J + JF is of course Lipschitz, 
with constant CA^-^+^^Z^. For the term FP + PF, we need to study P, given in (5.17): 

P = uj{p, kfQ + i(JF - FJ) - |(27r)=''^A2(ni(p, k) + n,{p, -k)). (11.22) 



We have 



uj{p, kfQ = T{2p) + (T * A){2p)uj{p, ky^ + pi + cu(p, kf 



r 



where pi collects the n > 2 terms coming from the expansion (8.16). From Proposition 9.2, 
the last term in (11.22) at J = r = is given by (9.12). Propositions 9.3, 9.5 and 9.6 
control the corrections. The (JF — FJ) term is bounded as in (11.21). To summarize, let 
us call Ps the sum of ^(JF — FJ) + uj{p,kyr, of the corrections to the J = r = term in 

— '^(2TiY"^\'^{ni{p,k) + ni{p,—k)^ and of the term corresponding to m in (9.12). Let Pe be 

the rest, i.e. T{2p) + (T * A)(2p)uj{p, k)^-^ + pi, and what corresponds to the sum in (9.12). 
Collecting the bounds established for these various terms, we get (remember that inverse powers 
of N are small compared to A^): 

Proposition 11.2. P can be written as 

P^Pe + Ps 

with 

Pe{p, k) = T{2p) + (T * A)(2p)a;(p, k)-^ + Pe{p, k) (11.23) 

where 

00 

Pe{p. k) = J2 Fni2p)hnip, k) (11.24) 

n=l 

where G E, with ||Pn|U < C'^II^IIe? f^f n>2, \\Fi\\e < CX^WAWe, and where the functions 
hn are smooth with 

WhnWoo < C", (11.25) 

Ps & S and 

\\Ps\\s < C{\\r\\s + X'{\\J\\s + \\t\\s + \\A\\e)). (11.26) 
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We may now bound P^(Pr + TP). The Ps term is in S, like J and can be bounded, as in 
(11.21), by the RHS of (11.26) times iV-i+5"/4. The T(2p) in (11.23) drops out from P^, since 
it is constant in k and we use (10.2). The other terms in (11.23) are bounded using 

ldq\F{q)\<C\\F\\E, (11.27) 

which follows easily from: 

I^WI S (Hp) + ^ + jviz^ + TvW^I^)"^"- ("-^^^ 

using (11.18) and: 

J \d{p)\-^dp < CN^-'^ k>l, (11.29) 

which is proven like (11.18). Thus P"'"(rPg; + PeT){p, k) is smooth, since, looking at (7.11) we 
see that the dependence on p, k in P-'-(rPE + PET){q, q + k — p) is only through the second 
argument of Pg, in which Pg is smooth, by Proposition 11.2. Moreover, using (11.27) and 
(11.23), (11.24), 

I j P^{TPE + PET){q,q + k-p)\<C-f\\A\\E. (11.30) 

Since P^(VPe + PE^){p,k) is smooth, let us identify it with an element of S of the form 
(0, 0,^, 0). Then, by (8.21), P^(rPE + PeT) e S and 

\\P^{rPE + PeHIIs < CN-^+''^^N'-''/^\\A\\e = CN~''/^\\A\\e. (11.31) 
Combining the above bounds, wc get: 

||P^(rP + Pr)||5 < CN-'^/^\\A\\e + C7V-i+W4Q|^||^ + ^py + ll^ll^) (1132) 
and P^(rPE + PeF) is Lipschitz in A, E, J, r, with constants 0{N-''/'^). 

We may summarize this discussion by rewriting equation (11.1); consider T,A given, for 
p & Eq and given for p ^ Eq. Let us denote 

sip) = d{-p){T * A){p) ^ d{-p)S{p). (11.33) 

Then 

( J„ = -V;^\l [(pit(2p) + p2s(2p), 0)^ + n] (11.34) 

for p e £^0) using the fact that = P'^Vg, and therefore, {Js,rs)^ = — H (Js,rs)^. Here, 
TZ e S includes all the terms in (11.1), apart from the first two in (11.9). Combining (11.12), 
(11.16), (11.21), (11.32), \\TZ\\s is bounded, for T, A, J, r as in the Theorem (using the fact 
that < 1), by: 

\\n\\s < CX'^BiT (11.35) 
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and is Lipschitz J, r with constant C(t), and in T, A with constant C(A^). ^ is the operator 
defined by (10.3), and bounded in Proposition 10.1. 
Let, for p e Eq, 

(Jo, rof = -V;'U(p,t{2p) + p2s{2p), Of (11.36) 
Then, by Proposition 10.1 and Lemma 9.1.b, 

||(Jo, ro)||5 < CX-'{\\t\\s + \\s\\s) < CX-\\\t\\s + \\A\\e). (11.37) 

Now, given Wg G B2 and Wi G B3 we have from Proposition 10.1. and (11.35) 

\\V-^Un\\s < CBiT, (11.38) 

i.e. V'^UTZ e B2 for B2 > CBi. It is Lipschitz in J, r, with constant C(t), and in T, A with 
constant O(A^). Thus we get, from the contraction mapping principle. 

Proposition 11.3. Given T,A,Wi as above, equation (11.34) has a unique solution {Js,fs) G 
B2, which is Lipschitz in (T, A) with Lipschitz constant C(A^) and in we with Lipschitz constant 
0{r) . 

This proves part (a) of Proposition 11.1. Now assume that we have proven part (b) of that 
Proposition. This will be done in the next subsection. We want to prove here part (c), since it 
is done in the same spirit as part (a). 

Thus, consider equation (11.3), for p ^ Eq. Following the argument that led from (11.1) to 
(11.34), we may rewrite it as: 

(J,r)^ = -V;^ [{p^t{2p) + p2s{2p),Qf + n\ (11.39) 

where Vp was defined in (10.4), for p ^ Eq, and is invertible by Proposition 10.1. We can write 
(11.39) as: 

( = -{Js, fsf - V;' [{p,t + P2S, + n] (11.40) 

where {Js,rs) is the Lipschitz function of w^, given by Proposition 11.3. The 5-norm of the 
RHS is bounded by CX~'^Bit and so, taking B^ — CBi, (11.40) is in B3, provided we show it 
vanishes for p E Eq. 

By assumption, T,A here is such that (11.2) holds for p E Eq (part (b) of the Proposition, 
to be proven below), and, by part (a), we know that (11.1) holds for p G Eq. Putting (11.1) 
and (11.2) together, we see that the full equation, (11.3), is satisfied by Qt,a + Wg + Wi for 
p G Eq. But (11.39) is merely a rewriting of (11.3), whenever Dp is invertible. Since we can 
assume, by choosing, if necessary, B in Proposition 10.1 larger than here, that "Pipo and X^Tripo 
are invertible, and since, by assumption, 

Qt,a + Ws + we^ Qt,a + Ws 

for p G Eq, we know that Qt,a + Wg solves (11.39) for p = ±po or p = tt ± po- But that means 
that the RHS of (11.40) vanishes for those values of p. We can define both sides to be zero for 
other values of p G E'o, if we want, without affecting the fact that J^, are in S, and thus we 
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obtain part (c) of Proposition 11.1, by applying the contraction mapping principle to the fixed 
point equation (11.40) to (Je^rc) e B3. 



Remark. The function J(T, ^4) is a general form of the Fourier Law, expressing the Eq^Py 
correlation function as a function of the local temperature and chemical potential. In particular, 
see (11.36) and remember that 11 is the identity of the first component, 

Mp, k) = Ki{p, k)t{2p) + K2{p, k)s{2p) (11.41) 

with 

Inserting to (7.7) and (7.14) we get 

{jo{p),Jo{p)f = ^^o{p){tip),s{p)Y (11.43) 

where k,o{p) is C*" in p. Since Dp is of order A^, ko{p) is of order A~^. The full K,{p) introduced 
in the Theorem has corrections 0{1) coming from applied to the parts of TZ in (11.34) that 
are linear in t, s. Since TZ is or order A^, the result is 0{1), i.e. small compared to kq, at least if 
the latter does not vanish. In the next subsection, we shall need the fact that K,{p) is invertible 
for p e Eq. To show that, let us first compute ko(0), which we shall denote for simplicity kP. 
Prom (10.3) and (6.8), we get 



£22(0) 



Inserting (7.5) to (11.10) and (11.11) we have 

Pj{0, k) = a sin ku!{k)~^ 

for J = 1,2, and so, defining 

'iljj{k) = 4sinA;a;(A;)~', j = 1,2, 

we have 

/^,(0,A;) = -z(Ai(0)-V,)(^)- 
Inserting to (7.7) we have, for j = 1,2, 

^-2 j dk sink uj{k){Cn{0)~%){k). (11.44) 

Prom (7.14) we get 

kI- = -2 j dk sink ri{Q,k)u;{k){Cu{Q)-^ipj){k). (11.45) 
where 7/(0, k) equals u;{k)''^ — J dkuj{k)~^. Now, note that 

< = -^(V'i,>Cr/(0)V',) (11-46) 
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where (•, •) is the scalar product in Tfo, and that 

4 = -^(^2,/:r/(o)^,) +/3o^(^i,/:n'(o)^,), (11.47) 

where /3o = / dkuj{k)~^. Computing the determinant of the 2x2 matrix k'^, we see that it 
equals the one with (3o = 0, and the latter does not vanish because £11 (0) is a strictly negative 
operator (see (10.6)). Now, from Holder continuity, we get that ko{p), and therefore K,{p), is 
invertible for p small, i.e. for p e £"0, for A small enough. 

Finally, observe that, since £11 (0) is a strictly negative operator, k{p), is a positive ma- 
trix, for p G Eq. To understand the connection with the Fourier law (1.1), note that t{p) = 
d{—p)T{p), and, in x-space, d{—p) is —V. 

11.2 Solving the conservation laws 

We are left with the proof of part (b) of Proposition 11.1. This reduces to solving the two 
conservation laws, equations (7.3) and (7.13), which are equivalent to (11.2). Indeed, (7.3) is 
(4.20) for X = y, which is the same as integrating (5.19) with dk, or taking the scalar product 
of (5.19) with uj{p,k)^^ in Hp. For (7.13), it amounts to taking the scalar product of (5.19) 
with a linear combination of a;(p, k)~'^ and uj{p, k)~^ . 

Let us introduce a more compact notation. We set J — and write (7.3) and (7.13) 

as 

d{p)J{p) + T{p) - e(p) = 27(Ti + e^^^Ta, 0)^ (11.48) 

for p E Eq with Eq given by (8.22), where we can assume that K{p) is invertible. Equation 
(11.48) has the friction term, see (7.15) 

^(p) = 7/ dkdqP(q/2,k)(l + e'^''+^^'')(2,^l;(p,k,q)f, (11.49) 

where we use e*^^ = e~*^-^, and the projection of N22' 

e{p)^{o,e{p)r, (11.50) 

where 9 is given by (7.9). 

J is given by the Fourier law i.e. the solution of (11.39) with leading term (11.43) which 
we may write as (see the Remark at the end of subsection 11.1): 

j{p) = K{pmp),s{p)f+z{p)], (11.51) 

where s is given in eq. (11.33). JT", and B are functions of T, A G E, and We G (including, 
indirectly, as functions of Ws, which is a function of T, A, Wi, by part a of Proposition 11.1). 
So, for wi fixed (11.48) is a nonlinear, nonlocal elliptic equation for T and A. We look for 
a solution to (11.48) in the ball Bi G E x E, defined in (11.4, 11.5). For such T the map 
A~^T*A = S is invertible (because T in (11.4) is close to a delta function, which is the 
identity for the convolution) and 

\\S-T+A\\E<Cr\ 
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Thus we can use T, S as the unknowns, in the ball Bi. The following Proposition collects the 
properties of the functions z and 9, studied in Propositions 11.3 and 9.7 (since |^(0)| < CN~^, 
we have that d-^{e - ^(0)) e S) : 

Proposition 11.4. z and d~\Q — 6(0)) are Lipschitz functions Bi x B3 ^ S ® S with 

\\z\\s < CX'^r, 

\\d-\e-e{0))\\s<CW (11.52) 
|e(0)|<CA^-^ (11.53) 

and Lipschitz constants bounded by CX^. 

Recall next that f E S is decomposed as (see (8.6)) 

f{p, k) = U{p, k)a42p) + f_{p, k)a42p), 
with cr±(p) = 1 ± e*^^, and f± = |/. Insert 

1 + e^(.+P)iv ^ 1 (^^^(p)^^(^) + a_(p)a_(g)) (11.54) 

into eq. (11.49) and use a+{q)a-{q) — ior q e ^Z2n to get: 

^{p) = ^+{P)(^+{P) + ^-{p)(^-{p), 

with 

J'±ip)^l j^dq J dkP±{q/2,k){2,'i/j{p,k,q)f (11.55) 

where J±dq = dqa±{q) i.e. the g-sum in /j_ runs over the odd (for — ) or even, non zero 
(for +), multiples of We used here the fact that c±(?)^ = 0,4 to cancel the factors of | in 
(11.54) and in P± = |P. 

Thus, (11.48) becomes two equations, one (+) valid on the even sublattice, and the other 
one (— ) on the odd sublattice: 

dj± + - e± = 27T±(1, 0)^ (11.56) 
with the average temperature (11.6) and 

T_ = |(Ti-r2) (11.57) 

It is useful to separate from (11.55) the part which is constant in p: 

^±(p)=^±(0)+/±(p), (11.58) 

and similarly for 0±(p): 

©±(p) = ©±(0) + e'±(p), (11.59) 
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Since if) in (7.16) is smooth we may write: 

f±ip)^d{p)g±{p). (11.60) 

and we may estimate g± as in the derivation of (11.31), and obtain (there is no P-^ here, so we 
have a term hnear in T) : 

Proposition 11.5 The functions g± & S are Lipschitz in {T,S) G Bi with ||5'±||s (ind the 
Lipschitz constants 0{N~°'^^). 

Using (11.51), we may write (11.56) as 

d{p)K{p) ((i±(p), s±{p)f + w±{p)) = 27r±(i, 0)^ - .F±(o) + e±(o) (11.61) 

where W±{p) = z±{p) + n{p)^^{g±{p) — d{p)^^Q'j_{p)), is defined for p G Eq, where K.{p) is 

invertible (see Remark at the end of subsection 11.1). So, combining Propositions 11.4 and 
11.5, W± G S, with 

||W^±||5 < CAV, (11.62) 

and the functions W± are Lipschitz in (T, S), with constants CX^. 

Let us next analyze ^±(0). Using Proposition 11.2, S = T * A and (7.8), we write 

^±(o)-e±(o) = 0± + V'± 

cp^^^l^dqj dk(T{q)+p{q/2,k)S{q)){2,,/;{0,k,q)f (11.63) 
ip± = 'yj^dqj dk(PE{q/2,k) + Ps{q/2,k)){2,ij{0,k,q)f-e±{0).{llM) 

For tjj^, we proceed as in the proof of (11.31). The contribution coming from Ps is smaU, 
see (11.21), while the one coming from Pe is to leading order, see (11.24), C(7A^t^). The 
contribution of ©±(0) is 0{N~^), by Proposition 9.7. So, V'i are bounded by 

|^±| < C7AV (11.65) 

and are Lipschitz in (T, 5*) with constant C'-fX'^. 

It is instructive to solve first the simplified equation (11.61) with W± dropped and ^±(0) 
replaced by (f)±. Then, for p 7^ 0, 

i^{p){t±{p),s±{p)) = d{p)-'C± (11-66) 

with 

e± = 27T±(l,O)^-0± (11.67) 
and, for p = (where obviously only the equation with index + holds), 

27T+(1, 0)^ = 0+. (11.68) 



49 



Equations (11.67) and (11.68) imply 

e+ = (11.69) 

so, by (11.66), 

(t+,s+) =0 

and, from (11.66) for — , we can write 

(t_(p),s_(p)) = rf(p)-^(r(p),C(p)) (11.70) 

with 

^{p){r{p),ap)f = C- (11-71) 

constant in p, i.e.: 

{r{p),ap)Y = K{p)-'K{0){T{0),aO)f, (11.72) 



where (t(0), C(0)j is defined by extending (11.71) (which was derived on the odd sublattice) to 
p — 0. Hence, since s{p) — d{—p)S{p), t{p) — d{—p)T{p), 

U{p) = {T{p),S{p)f = Uo6{p) + |ci(p)|-'(T(p),C(p)) V(p). (11.73) 

The unknowns are Uq = {Tq, Sq)^ , r(0) and ({0) and they will be determined from (11.67) 
and (11.68), where the functions given by (11.63), are functions of C/q, t(0), ({0) via 
(11.73): 

0+ = 7! dk{To + p{0,k)So) {2,ij{0,k,0)Y (11.74) 
= 27 J_dq J dk\d{q)\-'{r(q)+p{q/2,k)aq)) {2,^(0,k,q)f , (11.75) 

where the prefactor of 2 comes from the fact that, in (11.73), cr-ip) = 0, 2. 

Equations (7.8), (7.12) and (7.16) imply ij{0,k,0) = 2{p{0, k) - J dkp{0, k)), with p{0,k) = 
u{k)-\ So, 

Jij{0,k,Q)dk = 0, 
J p(0, k)dk = /3i > 0, 
I J p(0, k)i/j{0, k, 0)dk = J uj{k)-^dk - ( j ojiky^dkf = /32 > 0. 

Hence (11.74) gives 0+ = 27(To + /3iS'o, /?25'o)^ and, from (11.68), we get: 

^0 = , To^T+^Ut, + T2). (11.76) 
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To analyze (11.75), we need the straightforward 
Lemma 11.6 Let a be a function on j^1-2n- Then 

2 ljd{q)\-^a{q)dq = 7V(/±a(0) + 0(7V-")||a|U) 



;il.77) 



where /+ = ^ and /_ = 3/+ = | . 



The constants /+. /_ follow from the fact that (see (5.1)) 2 \d{q)\^^dq equals, to leading 
order in N, N times the sum over even or odd non-zero integers n (positive and negative) of 
^j^. The even sum equals a quarter of the sum over all integers, and the latter equals |. 



Inserting (11.72) into (11.75), and applying Lemma 11.6 to the result, we obtain 



2jNI_ 



/r(0) 
U(0) 



;ii.78) 



Recalling that = «;(0) (r(0), ((0))^ = 2^T_{l,0f - (see (11.71), (11.67)), we can write: 



2-fNI_ 



J A ) . 0(N-n] ( ) + .(0) ( ) ^ 2.r. ( J ) . (n.79) 



Since ^ — N i+'^Z^, we get: 



t(0) = (Ar/_)-i(T_ + C>(Ar-°/4)) 
C(0) = 0{N-^-'^'^). 



So, the simplified problem is solved by 



{T\ S') ^ (r+, 0)5(p) + ^|rf(p)|-V_(p) [(1, 0) + 0(7V-"/4) 



(11.80) 



In x-space, this is a linear profile: use, for p 7^ the explicit formula (8.4) with Jq replaced 
by = and observe that, in (8.3), j{x) = ^ for a; € [1,-^]; remember also that 

t{p) = d{—p)T{p), and, in x-space, d{—p) is —V. So, we get: 



T'ix)=T,+ ^-^iT,-T,) 



;ii.8i) 



plus a 0{N ^ "Z'^) correction to the slope 1/N. The first term comes from integrating (11.80) 
over p, using (11.77) for the second term, and Ti = T+ + r_. 

Let us now return to the full eq. (11.61) and incorporate the corrections W± in (11.61) and 
ip± given by (11.64). Let 



u±{p) = (t±{p),s±{p)y 
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and 



Write 



u±{p) = i^{p)-\d{p)-'C±) + v±{p). (11.82) 

In the absence of W±, 'il^±, v± — and ^± are as above. Consider the equation 

v±{p) = -W^ip). (11.83) 

Since the functions W± are Lipschitz in U, (11.83) has a unique solution v± & S, with 

ll^ills < C'AV, (11.84) 

and v± is Lipschitz in Uo,i±. With this v±, (11.61) becomes, as in (11.67) and (11.68) and 
(11.69): 



e+ = 0, 27T+(l,O)^ = 0+ + V'+, 
e- = 27r_(l,O)^-0_-V'- 



(11.85) 



Proceeding as in the simple case, 



where the term 0{'fTX'^) collects the contributions from v±- and that are bounded by (11.65), 
(11.84), and is Lipschitz in Tq,So,^±. Thus Tq,So,N^_ have C(rA^) corrections and 



\\{T,S)-{T',S')\\e<CtX' 
which, combined with (11.80), yields the claim of Proposition 11.1 (b). 



(11.86) 

□ 



A Holder regularity. 

In this Appendix, we prove some refinements of the Holder continuity of the kernels proven in 
Section 9. We start with a corollary of Proposition 9.3 that will be needed in Appendix B. For 
this, let G e C°(T^°') and consider the function 

g(p,k)^ I G(k)iip,k(dk) (A.l) 

on T^^*^, where jip^kidk) is defined by (9.43): 

3 3 

l^p,k{dk) = f y^ SiUj{ki) + S4U>{p — k) + iej ki — p — k)dkidk2dk^. (A. 2) 

1 1 
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Denote 



9 = m (A.3) 

Then we have 
Corollary A.l. 

(a)IfG IS smooth, I{G) is in ^"(1^+^). 

(h) Lei G{k) = H{ki, k2)h{kz) with H smooth and h e C°(T<^). Then I{G) is in ^"(1^+'^) and 

\\I{G)\U<C{H)\\h\U (A.4) 

(c) I is a bounded map from C°(T^'^) to C°(T^+'^). 

Furthermore, if G depends smoothly on some parameters so does g. 

Proof. By comparing (A. 2) and (9.31), we see that, if we replace k^ in (A. 2) by ks +p, we 
may identify ko, here and k' in (9.31). But then, the function G in (a) only changes in (9.31), 
and the statement (a) can be proven just as the one on M(p, k) in Proposition 9.3. In (b), the 
function H again affects only ps, and h is the function on which the operator Kp in Proposition 
9.3 acts. Hence, (b) follows from the claims made on Kp. Finally (c) follows because we can 
bound the integral I{G) by the sum norm of G and the resulting integral is continuous on T^^'^, 
for the same reason that M(p, k) is continuous. □ 

We need one more regularity result for the analysis of the 9 in (7.9) that will be made in 
the Appendix B. For this, define the function 

I{p) = J 6{u;{k,)+u;{k2)-u;{ks)-u;{k4))s(^p-Y,k^xik)dk. (A.5) 

Then 

Lemma A. 2. Let x in (A.5) be smooth and let I'{p) be the discrete derivative on ^ T-in- 
Then, for some a > 0, 

||/'||a < C (A.6) 

uniformly in N. 

For the proof we need a further Lemma: 

Lemma A.3. Let x(^) be smooth on and F{x) = sin(xi+X2+X3)— sinxi— sinx2 — sinxs. 
Then g{p) = J S{F{x) — p)x{x)dx is smooth if p and 

l5(p)|<c(logp)^ 

logp 



\9'ip)\<C 



p 
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Proof of Lemma A. 2. In (A. 5) write /cj = | + A;-, i = 1, 2 and = — | + k'^ which imply 
k4 = p — k[ — k'2 — k'^ — ^ and the argument of the 5-function equals 

3 

J2 sin(A;^ + k^ + k^) - sin k^ - sin k^ - sin k^ + P, 

a=2 



where 



Thus, 



with 



P = sin(ki + k2 + ka — p) — sin ki — sin k2 — sin k^. {■^■'^) 



/(p) = J dk j(P(k,p),k), (A.8) 



J(A,k) = J dxdyS(F{x) + F{y) + X)x{k,x,y) 
dt J dxg{t,k,x)5(^F{x)+t + Xy 



where 



g{t,k,x) = j 5{F{y) -tjx{k,x,y)dy, 



which by Lemma A. 3. is smooth in all variables, for t ^0, and is bounded by \g\ < C(logi)^. 
Similarily the x-integral yields 

J(A,k) = j dtds h{t,s,k)5{t + s + X), 

with h smooth ii t,s ^ and 

log 1^1 n I I 
(log|s| 



\dth\ < C 

and similarily for dgh. Then one gets 

|9,J| <C(log|A|)^ (A.9) 

Consider first the ^ 00 limit of (A.8). Then by (A.9) |/'(p)| < C/rfk| logP|^ 

Since P is an analytic function, | logP|^ is integrable. It is easy to do the argument for the 
Riemann sum. In the same vein, one can extract a little Holder continuity for I'{p). We leave 
the details for the reader. □ 

Proof of Lemma A. 3. a) We have 

daF — COs(xi + ^2 + X3) — COSXq,, 

so VF = <^=^ cosxi = COSX2 = COSX3 = cos(xi -\- X2 -\- x^) and thus xi — X2 — —x^ and 
permutations thereof. At these points F = 0. Hence F{x) 7^ =^ VF(x) 7^ 0. Thus, given 
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X e F ^(p) there is a neighbourhood U oi x and a smooth diffeomorphism: (pp : -Be(O) — > U 
such that (F o (pp){y) — p + yi and 0p is smooth in p. Given t/j e C^{U), 



9^ 



{p) = j S{F -p)'ijjdx^ J dy2...dyd{ipo(f)p^){0,y2,...,yd) detD(l)p 



is smooth. By a partition of unity, we conclude that g is smooth for p ^ 0. 
b) The Hessian of F is 

dadpF = Saf3 sin x^ - sm{xi + X2 + x^). 
At the critical points xi — X2 — —X3 it equals 

/Oil 

H^-sinxi 10 1 
V 1 1 2 

H has eigenvalues 0, sinxi, — 3sinxi. By a partition of unity argument it suffices to study two 
cases: 

1° X has support in a small ball around a critical point with xi 0, Xi ^ tt. 

2° supp X is a small ball around the origin, around (tt, tt, — tt) or permutations thereof. 

Case 1° By scaling and rotation, there exists a local coordinate z = {u,v,w) such that 

F ^uv + f{u,v,w) (A.IO) 

with / analytic, 0{z^) and /(O, 0, w) = = Du^„/(0, 0, w). Indeed, set Xi — x + u + w,X2 — 
x +v +w and xs — X — w, then 

F — sin x{cos{u + v + w)— cos{u + w) — cos{v + w) + cos w) + 
cos x{sm{u + v + w)— sm{u + v)— sm{v + w) + sin w) , 

which, upon scaling, is of the form (A.IO). Writing {u,v) = y, f{u,v,w) = {y,A{z)y) with 
A — 0{z). Since uv is nondegenerate, there exist an analytic diffeomorphism g close to identity 
such that F o g — uv. 
Hence 

9ip) = J 6{uv — p)x{u,v,w)dudvdw, 
with X e C^{B,{0)). Let ^ Jxdw, then, 

9{p) = - du\ogudu{(j){u,p/u) + (j){-u,-p/u)) 
Jo 

roo roo 

= -log|p| / dudu{'^{pu,l/u) - j du\ogudu{'^{pu,l/u)), 
Jo Jo 

with ip{u, v) = (p{u, v) + 4>{—u, —v). Thus g{p) = a{p) log \p\ + h{p) with a and b smooth. The 
claim follows. 

Case 2° We may suppose x = 0, the other points being similar. We have: 

F = S1C2C3 + C1S2C3 + C1C2S3 - S1S2S3 - si - S2 - S3 
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where Sj = sinxj, q = cosxj. The diffeomorphism near the origin = sinxj leads to 

F^{yi + y2){y2 + t/3)(yi + Vs) + 0{y^), 
and letting zi^yi + y2, Z2^y2 + ys, ^3 = ?/i + ?/3, we get: 

F = Z1Z2Z3 + f{z), 

with / analytic in B^{0), f — 0{z^), and symmetric under permutations of coordinates. We 
want to bound 

/ S{F{z)-p)tl;{z)dz, (A.ll) 
for ijj e C^(i?e(0)) and e small enough. By symmetry we may insert into (A.ll) 

6x(|^l| < \Z2\ < \Z3\)- 

Expand 

f{z) = fo{z2, Z3) + fi{z2, Z3)zi + f2{z)zl 

Since /i = 0{z2, z^Y, by a diffeomorphism close to identity we have (^2-2^3 + fi) ° <l> — Z2Z3 
i.e. may assume /i = and bound x from above by x(|2;i| < 2\z2\ < Sj^isl): 

< / S[ziZ2Z3 + foiz2,Z3)+f2{z)zl-p)'iP{z)-x{\zi\ < 2\z2\ < 3\z3\)dz, (A.12) 

with ijj e C^(i?e(0)). On the support of x and V-' I/2I < C'kal^ < C'e^ksl and thus | -22-23 + /2-2i| > 
(1 — (9 (e^) 1 2:2^3 1, since < 2|2;2|. So, given 2:2,2:3 and p, the argument of the S function in 
(A.12) either is nonzero for all zi e B^{0) or vanishes at zi{z2,Z3,p) satisfying 



P- fo 



Z2Z3 



< \zi\ < 2 



P- fo 



Z2Z3 



Moreover 



Z2ZZ + d^{f2zl)\ > (1 - 0{e^))\z2Z3\. Thus 



h{p)\<C I \z2Z3\-\{\z2\<2\zs\)x 



P- fo 



Z2Z3 



< Cej ijj{zi, Z2, Zs)dz2dz3. (A. 13) 



Since /o(-22,-23) = /o(0, 2:3) + 2:2/3 (-22, -23) where I/3I < Cz^ < Ce^\z3\, the second factor x in 
(A. 13) is bounded by 

P- /o(0,^3) 



X 



Z2Z3 



< Ce . 



Now, /o(0, Z3) is analytic, so 

/o(0,Z3) =<(l + 0(;^3)), 

for some a 7^ 0, n < cxd (actually n = 5). 

Insert this into (A. 13), write x = Xp{^2) + (x ~ Xp{^2)^, with 

Xp(^2) = x{\z2\ > C'IpI'^). 
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Then h = hi + h2 with 

hi{p)\ < J J^xiz2) J j^x(k3| > ^\z2\)xp{ze B, 
< C{\ogpf. 



For \z2\ < C|p|"n\ we note that, since, in the support of x ^ 

\p - a-23 I < Ce\z2Zs\, 
we may write Zs — (p/a)^/"" + x, with 

w-l 

" \x\ < Ce\z2\ IpI^/", 
where a is absorbed into C in the RHS, or 

II ^1 "-2 

\x\ < Ce\p\ " \z2\- 

So, doing the z^ integral, we get: 



P-/0(0,Z3) 
Z2Z3 



< 



Ce) 



dZ2 /, , ^ ^, ,_!/ 



/i2(p) <cej ^x(k2| < c^|p|^j|pr^/"k2| 

< Ce. 

We conclude that Isfd?)! < C(logp)^ i.e. the claim holds for g. The claim for g' is similar: dp 
brings an extra 1 2:22^3 and 



dz2dz'i / I ^ / n\ ^ 
^x(F2^3| > e{p)j < — 



[Z2Z3 



logp 



p 



□ 



Remcirk I(p) is not smooth for p 7^ 0. By some algebra one can show 

VkP = 0, P = 0^ki =k2 = k3, p = 2ki, 
and permutations. Thus zeros of P and Vj^P occur for nonzero p too. 

B Nonlinear estimates 



Here we prove the estimates on the nonlinear terms that were made in Section 10. They are 
based mostly on an analysis of convolutions of functions in E or 5", whose results were stated 
in Lemma 9.1. We start with the proof of this Lemma. 

Proof of Lemma 9.1 Recall first that j — j+o-+ + j+o-+ and T similarily. Let T * j — f. 
Then 
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Consider eg. T+ * and drop the +. 

Recall that j is defined through the 4-tuple j in (8.13). We need to define the 4-tuple j'. 
We set 

fo='^iT*m, (B.l) 
where the factor | follows from our conventions (8.13) and (6.5). Let, for p ^ 0, 

j[=T*n (B.2) 
j^ = T*j2 (B.3) 

T^7^ + T*^,-Ut*3i)- (B.4) 



(Ard)3/2 ^ {Ndf/"^ ' Nd Nd' 

Let us estimate these in turn, using two simple observations: 

II j fij)-p')gip')dp'\U<\\fU\9\W (B.5) 
where || — ||x,i is the norm, and 

II j fij^-p')d{p')-'dp'\U<\\f\U (B.6) 

which holds because the Hilbert transform of a Holder continuous function is Holder continuous 
(for a < 1, see [30], Theorem 106), and this is true also when is replaced by d{p)~^, and 
when we have discrete sums as here. 

Then, using ||r||ii < C||r||£;, (see (11.27)), and (B.5), we get: 

iij;iu<c^iiTiuiijiiu (B.7) 

lU^IU < ^^l|T||E|U2||a. (B.8) 

For jg, we use the identity 

d(p) = dip') + dip - p') + dip')dip - p'), (B.9) 

to write: 

(T*h-WdiT* Ji)) (P) - m i^diT * ^) - (T * ,0) ip) 

= Ml JiiP-p')id{p-pTW)T{p')dp' + 1 3i{p-p')d{p')T{p')dp') (B.IO) 

Hence, using (11.28), (B.IO) is bounded by 

C\\T\\e\\3UN\p\)-' /(^ + + ^^^)(1 + |p - p'\)-^)dp' (B.U) 

The integral is bounded by 

Clog|7Vp|(^ + Ar-iW2). 
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Combining with (B.ll) we arrive at the bound 

\iBAO)\<C\\T\\Msjj;^. (B.12) 
Similar calculations give the bound 



IT* ^-4^1 < IITIIe- ' 



ljVti|3/2 1 — ii^- |jV(i|3/2 ■ 
Combining with (B.12) we get 

lb3(p)lloo<c||r|U||i||5. (B.13) 

Finally (B.l) is bounded by 

I / \Tip)j{p)\dP<C\\T\\E\\j\\s (B.14) 

using (11.29), e.g. 

l{N'\p\Y'dp<C, 

Estimates (B.7) , (B.8), (B.13) and (B.14) give the claim (9.5). 
Let us next consider part (c). We set 

Jo = Uj*km (B.15) 



and, for p 7^ 0, 



j[ = 7Vd(^ ^h- + H,) + j^Uoh + hji) (B.16) 
j'^ = j^k2 + k* j2 - iV-^+"/2^2 * h (B.17) 
* + H2 + Ujoh + koh), (B.18) 



Js _ h 



where we write -^d * ^^^^3/2 + w * ^^^^^3/2 = Hi + H2, and the iJj's are defined after (B.23) below. 
To proceed, we need the following bounds: using (B.5), (B.6), we get that, 

hi(p) = j f{p- P')9{p')d{p')-^dp\ (B.19) 

satisfies 

||/ii||a<C||/|U||5|U. (B.20) 

This holds because we can write in (B.19) g{p') = g{p') — g{0) + fl'(O); since g is in C", \g{p') — 
g{0)\\d{p')-^\ < ||^||ab'r^+°, i.e. {g{p') - giO))d{p'y^ is integrable; so (B.20) follows from (B.6) 
applied to the g{0) term above and (B.5) to the gip') — g{0) term. 
Using (B.9), we then get that 

h^ip) = d{p) j f(p - p')d{p - p')-W)d(pr'dp' (B.21) 

59 



also satisfies 

\\h2\\a<C\\fU\gU. (B.22) 
Now, consider (B.16). The bound (B.22) applied to the first term implies that its C" norm is 

(^(^"')l|jl|UI|fcl||a. 

Now, use (B.9) to write 



^^(m * (ivd)3/2 ) - ii * + Nd * (Ndy/^ + S * (iv!)i/2 ■ l^^-^^) 

Let the sum of the first and third term in (B.23), plus the corresponding terms in Nd{-^* ^^^^^3/2 ) 
define NdHi and let H2 = {Nd)~^{j^ * (Nd)'^/^ ) P^^^ corresponding terms in * J^^^^aJ^- 
Since II^^^IIli < CN-^\\m^, we get from (B.5) that 



(^IU<0(iV-')||ji|UI|A;3||oo, 



which controls the contribution of the first term in (B.23) to NdHi. This bound holds also for 

the last term in (B.23), of course (which is even C(A^~^^^)||ji||a||^3||oo)- The terms :^(io^i+^oii) 
are trivially bounded, so we get: 

U\U<CN-'\\j\\s\\k\\s. (B.24) 
Using (B.5) and the fact that < A^~^ log A^||ji||oo for all terms in (B.17), we get 

U\\a<CN-'+^l\ (B.25) 

Consider now (B.18). Going back to the discrete sum (see (5.1)), it is easy to see that 
the first term is bounded by CA^~^|A'"p|~'^/^||j3||oo||/i;3||oo) i-e. that the contribution of this term 
to 1^3(^)1 is less than CA^'^lljllsH/cll^. Consider now {Ndy^{j^ * ^j2)^/2 ), contributing to H2. 

Going back to the discrete sum, it easy to bound it by CA^^^|A^]9|^^/^ log |A^p|. Finally, the 
bound on wijoks + kojs) is trivial and we get: 

\Mp)\ <CN-Hog{N\p\)\\j\\s\\k\\s. (B.26) 

Finally, using repeatedly (11.29), we get 

\fo\<CN-'\\j\\s\\k\\s. (B.27) 

Combining (B.24), (B.25), (B.26), (B.27), proves part c) of the Lemma, since the bound on 
IIj'II^ follows from the previous ones. 

Let us finally turn to (9.6). Note that, using (B.9), 

d{T*A)=dT*A + T*dA + dT*dA 

Since dT & S iiT & E the claim, for p 7^ 0, follows from (a) and (c). For p = 0, we simply 
apply (11.29) to all the terms in (T * A)(0). □ 

Using this Lemma, it is rather easy to give the proof of the main estimate of section 9. 

Proof of Proposition 9.4. (a) Apply (a) of Corollary A.l to (9.41) and (9.42) to get 

» 3 3 

^(P, k) = j gip, k,p) n Fi{2pi)S{2p -Y,'2pi)dp (B.28) 
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where g is smooth in p and in p,k. m may now be defined as an element of E as in Lemma 
9.1.b, the smooth function not affecting the bounds. 

(h) U m = u or a m — V and i 7^ 3 (in which case i — 1,2, and, by symmetry, we can choose 
i = 1) we have 

3 3 
^(P, k)= [fl F,{2p,){ [ G ki -pi)^ip,k{dkm2p -Y.2pi)dp 

with G smooth, and /i in <S. Using the representation (8.11) for /i we write G /i(pi, ki — pi) 
as a sum of terms, with singularities in pi; by Corollary A.l.b, we get that integrating each of 
those terms with Hp^kidk) gives rise to a function that is in p, k. Thus, 

j G fi{pi,ki-pi)fXp^k{dk) = f{pi,k;p,p2,P3) (B.29) 

where / is in 5 in the variables pi, k depending smootly on P2,P3 and in p. The convolutions 
with the Fj's can then be estimated as in Lemma 9.1. 
li m — V and i — 3 we need to study 

2 3 
/ n F,(2p,)/(P3,P - k-p3){[ Gnp,k{dk))5{2p - ^ 2p,)dp, 

i=\ •' 1 

where f & S. Shifting pa by p, this becomes 

2 3 
/ 9{p,k,p)llFi{2pi)5iJ2'^pi)Mp + ps, -k-ps)dp 



i=l 



here g, defined in (A.l) is, by Corollary A.l.a, C" in p, k and smooth in p. Performing the pi 
integrals for i = 1, 2, we get 

I F{p, k,ps)f{p + ip3, -k - \pz)dp^ (B.30) 

where F is C° in the first two arguments and belongs to £^ as a function of the third. We may 
proceed now as in Lemma 9.1. (a) to define and estimate the quadruple g in <S corresponding 
to (B.30). E.g. the component of index 1, see (B.2), is given by 

9i{p, k)^ j F{p, k,ps)fi{p + |p3, -k - ^P3)dps (B.31) 

where f\ is in both arguments. Again, since F is integrable in the third argument and in 
the others the integral is C° in p, k. The other components of g can be bounded as in Lemma 
9.1. a. 

(c) We have to specify again the quadruple m E S corresponding to m. Wc define mj = 
and m^ip, k) to be the integral (9.41) or (9.42) corresponding to m. We use Corollary A.l.c to 
do the ki integrals and estimate the pi integrals by brute force. Since only the singularity at 
Pi = (or, by periodicity at tt) matters, we need the easy bounds (remember that the variable 
p is discrete!) 

\Nd{p)\-^ * \Nd{p)\-^ < CN-^ \ogN\Nd{p)\-^''^ 

j^-i+a/2 ^ \Nd{p)\-^ < ClogA^A^-2+"/2 < CN-^+'^Z^ log N\Nd{p)\-^/^ 

^-l+a/2 ^ ^-l+a/2 < C N-^+''\N d{p)\-^''' 
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together with \Nd{p)\-^/'^ < \Nd{p)\-^ (since A^|p| > tt), which allows to use the first inequality 
here in order to bound the convolutions with the last term in (8.11), to conclude 



\m{p,k)\\^<C N-^+^M\El[\\fi\\s 

Ij^k 



which is the claim (9.46) since we defined m = (0, 0, 0, m). □ 
Finally, we prove the estimates on the function 9 defined in (7.9). 

Proof of Proposition 9.7 

By the (3 <-> 4) symmetry in (5.13), we get (leaving out the factor |(27r)^'^A^), 

s 1 

I - ^3) - I - ^4)] K,kidp dk) (B.32) 

and 

Kpk = (I] Sjo; {ki)+ie) ^ 5(2pi)5{p - 2^Pi)5{p ki)dpdk 

i it 

Then the [ ] in (B.32) equals 

2 f ^ - + (e^^ - 1) r{p, k) (B.33) 

with r smooth. The integral in (B.32) has singularities when 

J2siuj{ki) = 0. (B.34) 

Recall that (B.34) forces 

E«i = (B-35) 

Consider the s such that (B.35) holds in (B.32), and, replace [ ] in (B.32) by the first term of 
(B.33). We define 

^ s^y ft'^.. - ^) • ".^(^ «■ («-3^) 

By symmetry, we may replace S3 by ^ J2i=i and, by (B.35) also by — |s4. Again, by symmetry, 
Ssuj{k2) may be replaced by | Yli=i SiUj{ki). So, the parenthesis equals 2,J(kly^ Si=i ^i^iki), and 

the sum cancels the factor SiUj{ki) + ie^ in v' . Hence, (B.36) equals 

^i(p) = -|E li{Ws.^^,mP^)^^P-'^T.P^nP-T.h)dpdk (B.37) 
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We decompose 9 as 

^(p)=^i(p) + ^2(p)+P^3(p) (B.38) 

where 62 has the terms of (B.32) with X^Sj 7^ and the first term of (B.33) inserted, while ^3 
corresponds to the insertion of the second term of (B.33). 

Consider first di given by (B.37). Remember that 

Ws{Pi k — p) — Q{p, k — p) + isuj{k)~^J{p, k — p). 

The terms with an odd number of Q factors vanish by the s — > — s symmetry. Consider then 

the term hnear in J. We insert Q = Qo + ^ and start with the term with no r. After shifting 
the ki variables by 2pi (and using 2p4 = 0), we obtain a sum of terms of the form 

2 

/ T * A*"i(2pi)T * /l*^^(2p2) J(P3, A;3 +P3) ti^^ih + 2p,) + a;(A;,))"'""* • 
■' i=i 

uik^yMh + 2p3)-'5{2p^)S [p-2j2Pi)S (E ki) dk dp. (B.39) 

Now, use the fact that u!{ki + 2pi) — uj{ki) = (e^*^ — 1)0(1), which implies 

{u{h + 2p,)+u;{k,)y^"'' = {2u{k,)y'"'' +n,{e^'P' - 1)0{1). (B.40) 

For the second term on the RHS of (B.40), let us choose i = 1, which we can do by symmetry, 
and insert it in (B.39), to obtain, after integrating over ki, k2, k^: 

h{p) =mj (e^^fi - 1)T * A*^^{2p,)T * A*"^(2p2) J(P3, A;3 +P3) • 

/(p, k^)5{2p,)5 (p - 2 J2Pi) dpdks, (B.41) 

with / smooth. For the first term on the RHS of (B.40), we obtain: 

h{p) = J T* A*"i(2pi)r* A*"^(2p2)J(P3,A;3+P3) • 

f{p, ks)5{2p^)5 [p - 2 J2Pi) dpdks, (B.42) 

with / smooth and even in k^. Doing the /ca-integral, we get, for (B.41), 

/ Jip3,k3+P3)fip,ks)dk3 = g{2ps,2p) (B.43) 

where we used the 7r-periodicity of the result. (7 is in 5* as a function of ^3, depending smoothly 
onp, with \\g\\s < C'll-'^IU- For (B.42), we get, since J is odd in k^ and / even, that the integral 

vanishes if J{p3, k^ + P3) is replaced by J(p3, k^), and thus, since / is smooth, the integral can 
be written as: 

/ Jips,k3+P3)f{p,ks)dks^ J Jips,ks)f{p,ks-ps)dk3^{e^'P'-irgi2ps,2p), (B.44) 
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using again the 7r-periodicity of the result. ^ is in 5" as a function of ps, depending smoothly 
on p, with Hulls' < C|| JII5. We write, using the constraints S{2p4), 5 {p — 2X^pj), 

e2^P3 _ 1 = (e^P _ 1 + i)(e-2*^'i - 1 + l){e-^'P^ - 1 + 1) - 1. (B.45) 

Expanding the product we see that the integral (B.42) equals the sum of terms of the form 7i 
and of the form: 

h = d{p) j T * A*"i(2pi)r * A*"^(2p2)^(2p3, 2p)5(2p4)<5 {p -'^Y.P^ dp dh. (B.46) 

The integral in (B.46) is a convolution of two functions in E with one in 5, hence, by Lemma 
9.1, it is in S. The prefactor dip) cancels the d~^ , so that this contribution satisfies 

IM-'/2||5 < (C5„,+„,,o|M|5 + {C\\A\\eT^+^^)\\J\\s, (B.47) 

and also 72(0) = 0. 

Going back to /i, see (B.41), (B.43), we obtain: 

h{p) =nij f{pi)F{p2)g{p3,p)S{2p4)6 (p - Y,P^) dp (B.48) 

with f & S and F & E. So, we have a convolution of two elements of S and one of E, i.e. the 
convolution of two elements of 5". Going back to the definition (8.11), we see that we can write 

h^I[ + I'l + I'l' (B.49) 

corresponding to the i = 1, 2, 3 terms in the convolution of two elements of S that are 
bounded in part c) of Lemma 9.1. Prom that Lemma, we get: 

\n{p)\ < ]^^^^(^^<^n.+n„0||^||5 + (C||^||^)-+-)||J||5, (B.50) 

If we identify d{p)~^I[ with an element of S of the form (0,0,0,T*r), we get from (B.50) and 

{N\d{p)\)-y'<c, 

\\d-'l[\\s < iC6n,+n,44s + iC\\A\\Er^+^')\\J\\s- (B.51) 
Next, we get, also from Lemma 9.1.C, together with the definition (8.11): 

\\I"\\a < N-'+^{CSn,+n,ms + {C\\A\\Er+-')\\J\\s: (B.52) 

Now, identify d{p)~^I'( with an element of S of the form (0,T*r, 0, 0), we get, writing d{p)~^I'{ = 
{Nd(p))-^Nr^, that 

\\d-'ms < N-^+-{Cdn,+n„om\s + {C\\A\\Er+^')\\J\\s. (B.53) 

For I"', we use (9.10), and identify d{p)~^I'" with an element of S of the form (0, 0, 0,*). Since 

{N\d{p)\)-Hog{N\p\) < C, we get 

\\d-'l';'\\s < {C5n,+n„o\\t\\s + {C\\A\\Er^+-')\\J\\s. (B.54) 
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Combining these estimates, we get: 

IM-'/i||5 < {C5nM\\t\\s + {C\\A\\Er^+--)\\J\\s. (B.55) 
We also get, by (9.7): 

|/i(0)| < CA^-^ (B.56) 

The terms with one J and one r are sums of terms of the form: 

lT*A*-'{2p,)r{p2,k2-p2)J{P3,ks-ps)(uj{2p, - h) + uih))'^'''' 
u{h)-^u{k^)-H{2p^)5 [p-2j2pi)s{p-J2 k?j dp dk 
Doing the ki and k^ integrals, and shifting k2, k^, this equals 

h = Jt* A*-i(2pi)r(p2, k2)J{ps, ks)f{p, k2, ks,p)5(2p^)5 [p - J2Pi) dp dk2 dks (B.57) 
with / smooth. We have T * E E, r, J E S. Proceeding as with Ii, we get 

Wd-'hWs < {C\\A\\Er\\r\\s \\J\\s. (B.58) 

We also have, by (9.7): 

\hiO)\<CN-\ (B.59) 

In a similar way, we may analyse 62 i.e. (B.32) with X^s, 7^ and the first term in (B.33) 
inserted. We note that the terms that are odd in Q vanish since, for those terms, because of the 
s — — s symmetry, the measure is proportional to S(^SiU}{kij^, which vanishes for Z^Sj 7^ 0. 

Starting again with the term linear in J and with r = 0, it is given by 

I T * A*^'{2p,)T * A*^'{2p2)J{ps,ks)hip, + k,,k,)6{2p,)6 {p - 2J2p,) dpdk, (B.60) 

with h{p, k) = h{—p, —k) smooth. By oddncss of J in k, we may replace h by h{k^ + Ps, k-^) — 
h{kz — P3, ks) i.e., near p^ = 0, h is C(p3). Similarily, using J(p + tt. A; + tt) = J(p, k) and the 
27r-periodicity of h, 

J JiP3, k3)h{p3 + k3)dk3 ^^J ^(P3 - TT, ks) [h{k3+P3 - TT, /Cg + Tt) -/l(/C3 - (ps - Tt), /Cg + Tt)] 

i.e. (B.60) may be written as 

J T * A*"i(2pi)T * >l*"^(2p2)(e'^^^ - 1) J(p3, k3)hip3, k3)5{2p^)5 {p - 2 Y,Pi) dp dk3.{BM) 

Writing e^*^^ — 1 as in (B.45), and expanding the product, we see that the integral (B.61) equals 
the sum of terms of the form Ii and of the form I2, i.e. 

12 = d(p) I T * A*"^(2pi)r * A*^'(2p2)J(p3, k3)h(p3, k3)5{2p^)5 (p - 2 Y,Pi) dp dh (B.62) 
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with h smooth. The integral in (B.62) is in S and the prefactor d{p) cancels the d ^, so that I2 
has the same bound as in (B.47). Finally, the term with one J and one r is again of the form 
(B.57). 

The remaining terms in 9i and 62 are of type and Jr^. These are bounded by brute force 

by 

logA^M ,„„l„l|2 , II ,112 



mwrr+wjin (b.63) 



and considered as elements of S of the form (0, 0,0, Since \d\ ^ < CN we obtain by 
combining eqs. (B.55), (B.58), (B.62) and (B.63) 

+ ^2)||5 < cms + \\A\\e + ||r||5 + \\J\m\J\\s. (B.64) 

We still need to estimate p93{p) in (B.38), i.e. the contribution to (B.32) of the second term 
in (B.33), which we can write as: 

uJ'\h -P)- u;-\h) -(3^4)= pfih, h) + 0{p^), (B.65) 

with / odd. Consider the terms where Wg. = Qo, Vi. 
We get a sum of terms of the form 

J Y[Fn,{2pi)uj ^ '''{pi,ki-pi)s3S4Uj{k3)uj{k4) ^ 

[cj{k3 - p)~^ - coiks)-^ + cuik^ - p)~^ - uj{k4)~^] v\dp dk). (B.66) 

Write: 

a;-2-"'(p„ h - V^) = uj{k,)-^~''^ + (e^^^^ - l)n,0(l). 

Therefore (B.66) gives rise to two contributions: the one coming from uj{ki)~'^~"'^; after inserting 
(B.65) in the [-] in (B.66), and writing p = —id{p) + C(p^), this contribution can be written 
as, : 

^d{p){Fr,, * * F„3)(p)(/(p) - I{-p)) + 0{p'){Fn, * * F„3)(p), (B.67) 
with Fn = T* A*", 

I{p) = 1 5{u{ki) + u{k2) - ouiks) - uik4))S - E Ht)dk, (B.68) 

where is smooth, and / is odd (since /{k^, k^) above, and hence 0, is odd). By Lemma 9.1, 
{Fn^ * Fn2 * Fn.^){p) IS in E, so d{p){Fni * Fn2 * Fn^^^){p) is in S. By Lemma A. 2., the first term in 
the RHS of (B.67) multiplied by d~^ is in S and so is the last one, since C(p^)(F„^ *Fn2*Fn^){p) 
is N"^ times a C° function. The second contribution, coming from (e^*^* — l)niO(l), is of the 
form: 



3 

i=l 
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d{p)ni J n Fn,{2pi){e''''' - l)ij{p,k,py{dp dk), (B.69) 



where ip is smooth. The integral in (B.69) is in S with norm bounded by {C\\A\\e)^''\ The 
other terms in 6^ are simpler to bound, and we get: 

\\d-%\\s < C{\\t\\s+\\A\\E). (B.70) 

Of course, ^3(0) = 0. (B.64) and (B.70) together with (B.46), (B.56), (B.59), yield the claims. 
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